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SECTION  1 
INTRODUCTION 


Since  using  nuclear  devices  to  test  military  hardware  is  clearly  out  of  the  question,  an 
alternate  approach  is  needed.  To  this  end,  plasma  radiation  sources  were  designed  and 
built,  and  one  such  device  is  the  Phoenix  plasma  radiation  source  (PRS). 

The  goal  of  the  Phoenix  project  is  to  produce  x  rays  with  wavelengths  less  than  1  A  to 
best  simulate  the  radiation  produced  by  nuclear  devices.  The  Phoenix  machine  produces 
the  radiation  by  exploding  wires  and  driving  them  into  a  z-pinch.  During  the 
compression  of  the  plasma,  electrons  are  stripped  from  the  atoms,  resulting  in  highly 
ionized  atoms.  To  produce  high  energy  x  rays  it  is  necessary  to  remove  all  but  a  few 
electrons  from  the  atoms.  Since  the  remaining  electrons  will  be  tightly  bound  to  the 
nucleus,  the  radiation  they  produce  when  they  undergo  transitions  will  be  in  the  x-ray 
region.  The  energy  of  the  x  rays  scale  with  atomic  number  Z,  and  it  is  necessary  to  use 
iron  or  copper  to  achieve  radiation  with  wavelengths  less  than  1  A.  If  copper  or  iron 
wires  are  used  in  the  Phoenix  device  however,  the  highest  ionization  states  of  the  atoms 
are  not  reached.  An  instability  is  thought  to  develop  which  prevents  a  total  collapse  of 
the  plasma  column.  There  are  several  possible  instabilities  which  can  develop,  but  it  is 
still  unclear  which  one  is  responsible  for  the  disruption  of  the  pinch. 

One  theory  is  that  this  instability  is  due  to  anomalous  electric  fields  (AEFs).(1)  These 
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fields  are  manifested  as  waves  which  are  disruptive  to  a  smooth  collapse  of  the  plasma. 
AEFs  have  been  experimentally  verified  in  many  plasma  devices  with  lower  densities 
than  Phoenix.(2)  It  is  suspected  that  these  AEFs  may  also  occur  in  a  dense  plasma  like  in 
the  Phoenix  device,  and  to  test  this  hypothesis  a  method  for  measuring  such  electric  fields 
has  been  developed.  The  effect  of  these  electric  fields  on  emitting  atoms  is  quite 
complex.  By  knowing  how  these  electric  fields  broaden  the  spectral  lines  produced  by 
the  atoms  in  the  plasma,  a  measure  of  the  electric  fields  can  be  made.  This  dissertation 
describes  and  discusses  the  experiments  performed  on  the  Phoenix  device  to  measure  the 
levels  and  angular  distribution  of  AEFs  in  a  superdense  plasma.  An  aluminum  plasma 
was  used  for  the  study,  since  the  hydrogen-like  state  can  be  reached  where  the  affects  of 
oscillating  electric  fields  on  line  broadening  are  well  known.  The  presence  of  these  fields 
in  a  plasma  of  this  density  would  suggest  that  this  type  of  instability  may  be  present  in 
superdense  plasmas  and  needs  to  be  considered  as  a  possible  explanation  for  the  inability 
to  reach  the  highest  ionization  states  in  higher  Z  elements. 
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SECTION  2 

DENSE  PLASMA  MODEL 

2.1  PRELIMINARY  CONSIDERATIONS. 

The  key  to  plasma  spectroscopy  is  understanding  the  internal  state  of  a  plasma  by 
observing  the  radiation  emitted  from  the  plasma.  To  do  this  it  is  necessary  to  have  a 
model  which  can  predict  lineshapes  based  on  the  physical  parameters  of  the  plasma.  The 
model  used  in  this  analysis  assumes  that  the  observed  lineshapes  result  from  three 
primary  broadening  mechanisms:  Doppler  broadening,  Stark  broadening  and  self 
absorption.  Zeeman  splitting  is  also  considered  as  a  possible  broadening  mechanism. 

In  the  Phoenix  device  it  is  expected  that  the  plasma  has  an  electron  density  of  Ne » 

1 0  cm'  ,  and  an  electron  temperature  of  Te  a  T;  »  1  keV.  For  this  density  and 
temperature,  the  electron  collision  frequency  coe  will  be  on  the  order  of  1013  s'1.  Waves 
that  do  not  have  a  frequency  higher  than  this  collision  frequency  will  not  propagate.  The 
only  waves  which  can  have  a  frequency  higher  than  this  collision  frequency  will  be 
Langmuir  waves  with  frequencies  at  cope  «  10b  s'l.(3)  For  that  reason,  the  source  of  the 
oscillating  electric  fields  in  the  model  was  taken  to  be  Langmuir  waves. 

2.2  LINE  BROADENING  MECHANISMS  IN  THE  PHOENIX  PLASMA  SOURCE. 
This  section  outlines  the  model  used  in  the  analysis  of  the  spectra  obtained  on  the 
Phoenix  device.  Reviews  of  the  theory  behind  the  broadening  mechanisms,  relevant  to 
the  model,  are  presented  here. 
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2.2.1  Doppler  Broadening. 

Doppler  broadening  of  spectral  lines  results  from  the  relative  motion  of  emitting  atoms  or 
ions  towards  or  away  from  the  observer.(4)  Photons  which  are  emitted  with  a  fixed 
wavelength,  as  given  by  the  difference  in  energy  levels  of  the  atom,  are  shifted  when 
observed  in  the  rest  frame  of  the  laboratory.  The  two  sources  of  ion  motion  in  plasmas 
are  thermal  motion  and  bulk  plasma  motion. 

For  the  thermal  motion  of  the  ions  where  the  velocity  distribution  is  assumed  to  be 
Maxwellian,  the  number  of  atoms  dN  with  velocity  between  v  and  v  +  dv  will  be  given 

by(5) 

dN(v)  =  N0  exp(-  •  (2.1) 


Using  the  Doppler  formula  X  =  X0  [1  ±  (v  /  c)] ,  the  expression  for  the  Doppler  broadening, 


which  gives  the  number  of  atoms  emitting  with  wavelength  between  X  and  X  +  dX ,  is 


dN(X)  = 


Me 2 


2t iKTX 


exp| 


o' 


me 


2KT 


\X0 ,/ 


dX, 


(2.2) 


Since  the  intensity  is  proportional  to  the  number  ot  emitters  in  the  range  X  to  X  +  dX , 

Eq.  (2.2)  gives  the  Doppler  broadened  profile.  The  resulting  half  width  at  half  maximum 
of  this  profile  is 


AX 


HWHM  ~ 


X  \2KT\og2 
2c  L  M 


(2.3) 
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For  an  estimated  ion  plasma  temperature  in  the  Phoenix  device  of  T»  1  keV  at 
X  =  6.053  A  (A1  XIII),  the  resulting  HWHM  is  1.35  mA. 

In  many  laboratory  plasma  devices,  such  as  tokamaks  and  stellarators,  the  bulk  motion  of 
the  ions  is  small  compared  to  the  thermal  motion/6,75  In  z-pinch  devices,  however,  the 
bulk  plasma  velocities  are  on  the  order  of  0. 1  %  c  and  are  generally  larger  than  the  ion 
thermal  velocities  due  to  the  presence  of  strong  driving  currents/85  The  bulk  plasma 
motion  then  plays  a  dominant  role  in  the  Doppler  broadening  of  the  spectral  lines.  Our 
model  assumes  that  the  velocity  distribution  of  the  bulk  plasma  is  Gaussian,  as  commonly 
done  when  there  is  no  information  about  the  velocity  distribution. 

The  velocity  distribution  for  the  bulk  plasma  motion  is  constructed  by  replacing  the 
thermal  velocity  by  an  average  bulk  plasma  velocity  v  such  that 


2 KT 
M 


(2.4) 


Introducing  a  new  parameter  b  =  103  v/c,  the  width  of  the  Doppler  profile  AXD  is  given 


by  AXD  =  10  3bX0(mA).  The  normalized  line  profile  then  becomes 


_  rAL 

G( AX,A\D)  =  -!-re  ,  (2.5) 

AXDy/n 

with  a  halfwidth  at  half  maximum  given  by  AX  HivHiw{m  A)  =  &X0(mA)lO~3  VIn2  .  A  b 

value  of  1  for  X0  =  6.053  A  leads  to  a  HWHM  ~  5.0  m  A.  This  value  of  b  corresponds  to 
a  bulk  velocity  of  0.1%  c. 
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The  complete  Doppler  broadening  of  the  line  will  realistically  be  a  convolution  of  the 


Doppler  broadening  from  both  the  thermal  motion  and  bulk  motion.  Since  both  profiles 
are  Gaussian,  the  convoluted  line  will  also  be  Gaussian.  The  width  of  the  convoluted 
profile  will  be  given  by  the  square  root  of  the  square  of  the  widths  of  the  individual 
Gaussians.  Since  the  broadening  by  bulk  motion  is  significantly  larger  than  that  of  the 
thermal  broadening,  the  resulting  width  will  be  close  to  that  of  the  bulk  plasma  motion 
alone.  For  this  reason,  the  thermal  Doppler  broadening  is  ignored  in  the  model. 

2.2.2  Stark  Broadening. 

Stark  broadening  results  from  the  change  in  energy  levels  and  wave  functions  of  an  atom 
due  to  the  presence  of  electric  fields.  In  a  plasma,  static  macroscopic  electric  fields  are 
not  present  due  to  the  shielding  effects  of  the  plasma.  The  atoms  encounter  electric  fields 
only  from  neighboring  charged  particles  and  from  waves  traveling  in  the  plasma.  We 
start  our  discussion  of  Stark  broadening  by  briefly  summarizing  the  solution  of  the  linear 
Stark  effect  using  parabolic  coordinates/9"1 1) 

2.2.2.1  Linear  Stark  Effect.  The  change  in  energy  levels  of  an  atom  due  to  a  constant 
external  electric  field  can  be  found  from  a  solution  of  the  Schrodinger  equation  as 
H  =  Hq  +  H' ,  where  H0  is  the  hydrogenic  Hamiltonian  and  H'  =  qFZ  is  the  interaction 
term  with  the  electric  field  F.  The  problem  becomes  one  of  solving  Schrodinger’s 
equation  for  the  system 


V2u  + 


2m 

fi 


Ze1 

E  +  +  Fez 


u  =  0 . 


(2.6) 
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As  long  as  the  electric  field  is  not  too  strong,  the  potential  term  arising  from  F  can  be 
treated  with  standard  perturbation  theory  using  the  wave  functions  for  the  unperturbed 
hydrogen  atom.  Schrodingers’  equation  for  a  H-like  atom  without  an  electric  field, 
solved  in  parabolic  coordinates,  has  the  well  known  solutions 


where 

e  =  ^E 

L £  =  Laguerre  polynomials 

n,  nh  n2,  m  =  parabolic  quantum  numbers  that  must  satisfy  n  =  n,+  n2  +  m  +  1 
n  =  principle  quantum  number 

m  =  magnetic  quantum  number  having  positive  values  less  than  n 
With  the  electric  field  present,  a  perturbation  is  performed  and  the  corrected  energies  to 
first  order  are 


E 


n,n{,n2 


1  Z2  Me4  3  ea0Fn 
2~nT~nT~  +  2  Zh 


(«i  ~n2). 


(2.8) 


The  second  term  is  the  linear  Stark  shift  of  the  energy  level  given  by  the  quantum 
numbers  n,  nb  and  n2. 


2.2.1  Allowed  Transitions  for  Lyman  |3,  y,  5  and  s. 

For  transitions  between  levels  in  hydrogen-like  atoms  the  dipole  selection  rules  must  be 
satisfied,  that  is  Am  =  0  for  the  n  polarization  and  Am  =  ±1  for  the  a  polarization/ 12)  The 
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n  polarization  corresponds  to  the  electric  field  of  the  radiation  along  the  direction  of  the 
applied  electric  field  and  the  a  polarization  corresponds  to  the  electric  field  of  the 
radiation  perpendicular  to  the  direction  of  the  applied  electric  field.  In  this  experiment 
we  are  interested  in  the  Lyman  y,  5,  and  £  lines  corresponding  to  transitions  to  the 
ground  state  (n  =  1 ).  Figure  2-1  shows  the  Stark  components  for  each  of  the  allowed 
transitions  with  corresponding  intensities  and  polarizations  for  each  of  these  lines.  The 
energy  shift  of  each  component  is  given  by 


2  Zh 


where 


Ak  =  n'(nl  -n2  )-n(n]  -n2). 


2.2.2  Electric  Fields  Arising  from  Charged  Particles  in  a  Plasma. 

The  interactions  of  ions  with  the  neighboring  charged  particles  depend  upon  the 
properties  of  the  perturbing  p articles. (4'13)  For  slowly  moving  particles  such  as  ions,  the 
electric  fields  do  not  vary  rapidly  in  time.  The  resulting  electric  field  can  then  be 
calculated  by  taking  a  statistical  average  over  the  perturbers  using  a  quasi-static 
approximation.  Fast  moving  particles,  such  as  electrons,  interact  with  the  atom  by 
collisions.  When  electrons  collide  with  the  ions  they  cause  phase  and  amplitude 
modulations  of  the  radiating  oscillator  and  shorten  its  lifetime.  All  three  of  these 
electronic  effects  lead  to  spectral  broadening  and  are  treated  with  an  impact 
approximation. 
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The  range  of  validity  for  these  two  approximations  depends  on  the  number  of  particles  in 
the  Weisskopf  sphere  given  by  g  =  (pw  /  p  v)3  where  the  Weisskopf  radius  pw  and  the 


average  distance  between  particles  p#  are (14,  ,5) 

3  nane 2 
P"  =  2  — 

and 


(2.10) 


Pw  =  iV-1/3.  (2.11) 

The  Weisskopf  radius  gives  the  effective  optical  cross  section  for  the  collisions  of 
particles  with  the  radiating  ion.  If  the  number  of  particles  in  the  Weisskopf  sphere  is 
small  so  that  g  «  1 ,  then  the  nearest  perturbing  particle  causes  the  dominant  broadening 
and  the  impact  approximation  can  be  used.  If  the  number  of  particles  is  large  so  that  g  » 
1,  then  the  atom  interacts  with  a  large  number  of  particles  allowing  for  a  statistical 
approach  and  the  quasi-static  approximation  can  be  used. 

2.2.2.1  Ion  Broadening  by  the  Quasi-Static  Approximation.  We  first  consider  the 
emitter  interacting  with  neighboring  ions.  Since  the  ions  are  massive  and  slow,  the  ions 
can  be  treated  as  fixed  for  the  calculation  of  the  electric  fields  they  produce.1  i6) 

We  need  to  calculate  the  distribution  W0(F)  of  electric  fields  F  that  an  emitting  atom  or 
ion  encounters  from  neighboring  ions.  We  assume  that  the  ions  are  statistically 
independent  so  that  the  probability  of  finding  a  field  F  in  the  interval  from  F  to  F  + 
dF  is  given  by 
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(2.12) 


f  n  ' 


Wn 


(F)  =  f"l 5  F~HFj  P{rx,r2,...,rn)dlrxd*r2,...drn2 


V  y=l  y 


where 


Z  er 
/r  --  £  i 


rj 


(2.13) 


the  Coulomb  field  of  each  neighboring  ion.  P(rx,r2,...,rn)  is  the  probability  of  finding 
ion  1  in  the  range  rx  to  rx  +  d  F, ,  finding  ion  2  in  the  range  r2  to  r2  +  d  r2,  etc.  Since  the 
particles  are  statistically  independent,  this  probability  will  be  given  simply  by  V~n. 

For  isotropic  plasmas  Wq  (F)  only  depends  upon  the  magnitude  of  F ,  so  we  have 

W(F)  =  4nF2  W0(F) .  Carrying  out  the  integration  we  arrive  at  the  result  first  obtained  by 
Holtsmark,  namely, 


W(u)  -  —  [  xe  x  2  sin(ux)dx , 
71  o 

where 

u  =  F/F0 

F0  =2.6  Z',3eNe2l3.(il) 


(2.14) 


For  our  model,  we  define  a  dimensionless  variable  C  for  density,  such  that  Ne=  C  1.56  x 
1022cm'3. 

To  get  the  line  profile  from  this  distribution  of  electric  fields  we  find  the  wavelength 
change  AX  associated  with  the  Stark  shift  due  to  the  electric  field  F, 


AX  = 


3  ccinXft 
4n  hZc 


n(n j  ~n2)F. 


(2.15) 
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For  a  transition  between  two  levels,  where  both  levels  are  split  under  the  electric  field,  the 
change  in  wavelength  from  the  unperturbed  transition  will  then  be  given  by 


where 

Ak  =  Stark  constant  associated  with  the  shift  between  the  two  energy  levels  given 
by  Ak  =  -n2)-n'(nl  -n2). 

Defining  A0  =  3ea0X2o  /  AnhcZ ,  then  F  =  A Xn^n,  /  A0  Ak  and  u  =  AXn^n,  /  A 0F0AJc,  and 
the  line  profile  for  ion  broadening  is  given  by 

/( AX,Ak)  =  jdu  W(u) 5(AX  -  A 0F0uAk)  = - - - W( AX  /  A0  F0  Ak)  .  (2.17) 

A0  F0  Aa: 

This  profile  does  not  exhibit  maximum  intensity  at  the  line  center,  see  Figure  2-2.  This 
occurs  because  the  most  probable  electric  field  that  the  atom  will  encounter  will  be 
different  from  zero,  thus  leading  to  a  shift  of  wavelengths  from  the  zero  field  values. 
Figure  2-2  gives  the  shift  and  the  width  of  the  Holtsmark  profiles  as  a  function  of  C  Ak 
for  the  Lp  and  L,,  lines  of  A1  XIII.  From  this  graph  we  can  get  an  idea  of  the  size  of  the 
broadening  due  to  the  neighboring  ions  as  a  function  of  the  density  and  the  Stark 


constants. 

2.2.2.2  Electron  Broadening  by  the  Impact  Approximation.  For  the  impact 
approximation,  we  require  that  the  duration  of  the  collision  between  the  electron  and 
atom  is  small  compared  to  the  time  between  collisions.  The  assumptions  of  this  theory 
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Wavelength(mA) 

(a)  Holtsmark  profile  for  Lyman  y 
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(b)  Wing  center  and  half  height  of  the  Holtsmark  profiles 


2/3 

C  Ak. 


Figure  2-2.  Holtsmark  profiles  as  a  function  of  C 


Ak. 


(a)  Relative  motion  of  the  atom  and  perturbing  particle  is  quasi-classical,  so  that 
one  can  use  the  concept  of  trajectories. 

(b)  The  trajectory  is  rectilinear. 

(c)  Interactions  with  the  nearest  perturbing  particle  play  the  principle  role  in  the 
broadening. 

With  these  assumptions,  the  resulting  profile  of  an  isolated  line  in  the  impact 
approximation  is  Lorentzian  and  the  profile  is  given  by 


j _ i 

(co  -  A)2  +(^)2 


(2.18) 


where 

y  =  2  JV'(vct')  ,  the  width  of  the  profile 
A  =  N(va”) ,  the  shift  of  the  profile  (19,20) 
a' ,  ct"  =  cross  sections  associated  with  the  width  and  shift. 

Physically  the  broadening  is  predominantly  due  to  the  disruption  of  coherence  of  the 

emitted  radiation,  so  that  the  radiation  from  the  atom  is  broken  up  into  wave  packets  each 

with  internal  frequency  co0.  The  Fourier  analysis  of  the  wave  packets  produced  leads  to 

broadening.  Under  the  impact  approximation,  the  width  of  the  resulting  line  y  is  given  by 

the  electron  impact  width 

7  e=JVeVeawlnA,  (2-19) 


where 


aw  =  Tip/ ,  the  effective  optical  cross  section  of  the  atom 
The  In  A  term  comes  from  the  limits  of  the  impact  parameter  in  determining  the  cross 
section/- n  The  electron  impact  width  in  terms  of  plasma  parameters  is  given  by 
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Y 


(2.20) 


vh.  ^ 


e  rjJ  2  £  2 

Z  ?n„  v 


e  K7e 


n  hu) 


pe 


where 


Te 


m,  . 


For  A1  (Z  =  1 3)  with  an  electron  temperature  of  1  keV, 

ye  =  3.7xl012«V.  (2.21) 

The  resulting  HWHM  of  the  broadening  for  the  Lyman  P,  y  and  5  lines  corresponding  to 
n  =  3.  4  and  5,  are  0.3,  0.8,  and  2.0  mA  respectively.  Since  this  broadening  is  small 
compared  to  the  other  broadening  mechanisms,  we  can  ignore  impact  broadening  by 
electrons  in  the  model. 

2.2.3  Oscillating  Electric  Fields  (Waves). 

We  next  consider  the  effect  that  oscillating  electric  fields  have  on  an  emitting  ion.(22)  The 
Hamiltonian  for  an  atom  in  presence  of  an  oscillatory  electric  field  can  be  written 

H  =  H0+H',  (2.22) 

where 

H0  =  the  Hamiltonian  of  the  unperturbed  H-like  atom 

H'  =  ZeEo  cos(cof),  the  contribution  from  a  monochromatic  oscillatory  electric  field. 
From  Schrodinger’s  equation 

¥  —  =  (H0  +  eZE0  cos(oj/))vF  .  (2.23) 

dt 

We  look  for  a  solution  of  the  form 


A  = 

n,k 
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where 


e~yt=  the  time  dependent  coefficient 

bnk  =  Fourier  coefficients 

\\f„~  the  hydrogenic  wave  functions. 

Inserting  this  expression  into  Eq.  (2.23)  and  using  the  orthogonality  of  the  hydrogenic 
wave  functions  we  obtain 

(r  ~En-  k(£>)bnk  =  Vnn  {bnJc_x  +  bn  k+l ) ,  (2.24) 

where 

Vnn  =Wn\^~\Vn)-  (2-25) 

We  consider  two  cases  for  the  solution  of  this  equation.  First,  if  y  is  different  from  En, 
then  (  y  -E„)  will  be  large  compared  to  the  perturbation  potentials  Vnn  and  bnk  =  0.  When 
y  =  En,  the  term  on  the  left  of  Eq.  (2.24)  will  no  longer  be  small  compared  to  the 
perturbation  potentials  and 

—kobnk  —  Vnn{bnk__  |  +bn,k+ jj.  (2.26) 

The  coefficients  bnk  then  satisfy  the  same  recursion  relations  as  for  Bessel  functions, 
namely, 

2k 

—Jk{x)  =  Jk+\{x)  +  Jk-\{x),  (2.27) 

x 

so  the  solution  to  the  problem  is  given  by 

Y  „  =  En ,  and  bnk  =  Jk  f— 1 ,  (2.28) 

v  co  / 

16 


with  the  wave  functions  becoming 


Wn=e~^kfe^Jk(^n. 

*=-CO  V  ©  ' 


(2.29) 


The  profile  of  the  line  is  calculated  from  the  wave  function  by 


Km  « 


2  JIT 


I  exp(-z7Aco)(vF/!  (x,  t)  j<x)|  x¥m  (x,  t))dt 


(2.30) 


where  d(x)  =  er  is  the  dipole  operator.'5,23)  Working  out  the  dipole  expression  first,  we 
have 

k~x  (  ?  v  ^ 

('?.  wi  'p.  > = *'<£-‘£-1'  5>,  kw  - 


J=X 


>=- “> 


2V_. 


j\  ~ '.m. 


CO 


(2.31) 


Inserting  this  expression  into  the  integral  given  by  Eq.  (2.30)  and  integrating  over  all  time 


lim  —  J  e-‘>[*<»-iE-Em+(j-k) CO)]  *y  “y  r  f  2 V„n)  T  (2Vm 


t->°o  2nx 


k=- oc  y=— oo 


X  .2 


=i"f  <2-32) 

1  £=-co  y=— oo  10 


The  line  profile  is  then  given  by 


k= °o 

■z 

k=—a 0 


5[Aco -(£„  -Em  +(/-£)©)] .  (2.33) 


Making  the  substitution  p  =  £  -j ,  the  profile  can  be  rewritten 


p=-ooly=-oo 

Using  the  following  relationship  for  Bessel  functions 

k= oo 

Jn(x~y)=  YJJk+n(x)Jk(y)’t or  n  =  0,  ±1,  ±2, 

£=-00 


S[Am-(£,  -£. -p<o)].  (2.34) 


(2.35) 


P=CC 


the  integral  simplifies  to 

I,im  (A®)  °c  |(v  „  |^(x)|  y  m  )|  ^  Z 

where 


J, 


(2Vim 

2V  ' 

mm 

l  © 

© 

5[Acd  - (En  - Em  -/>©)] ,  (2.36) 


Kn  =  (V„  |^L[V»)  =  ^^^-«(«l  -«2)' 


^  i  i  /i  /  «  vi  x  ✓  (2.37) 

2  2  2Zeme 

Taking  the  relative  frequency  change  to  -  En ,  we  have  the  expression  for  broadening 
by  monochromatic  electric  fields, 


L,(Ao))  |{v„  |^(x)[v|/m)|2  fVj(A^£0)5[A©  +/*»] ,  (2.38) 

P=—CC 

3h 

where  we  have  defined  £  = - .  This  calculation  shows  that  under  the  action  of  an 

2<S)Zeme 

oscillating  field,  a  line  is  split  into  a  large  number  of  satellite  lines  located  at  A©  =  pay, 
with  the  amplitude  of  each  satellite  given  by  J2p  (M^E0) .  We  refer  to  \p\  as  the  satellite 
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number,  since  p  =  0  corresponds  the  center  of  the  line,/?  =  ±1  to  the  first  satellite  located 
at±co,  etc. 


2.2.3.1  Superposition  of  Electric  Fields.  In  the  plasma,  radiating  ions  will  encounter  a 
large  number  of  waves  with  random  phase.  The  resulting  electric  field  E(t)  encountered 
by  the  ions  will  be  the  sum  of  all  of  these  waves,  that  is, 

N 

E(t)  =  Y,  EJ  cos(co t  +  (p; )  •  (2.39) 

j=\ 

We  can  rewrite  E(t)  in  terms  of  an  effective  electric  field  E0  and  phase  factor  (p  as 

E{t)  =  E0  cos(co t  +  (p) ,  (2.40) 


(2.43) 
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where 


1/2 


E0  = 


HEj 


(2.44) 


v  j  J 

In  most  plasmas,  Langmuir  oscillations  develop  anisotropically.(2)  With  the  axial 
symmetry  of  the  device,  it  is  expected  that  the  value  of  the  electric  field  parallel  to  the 

current  E0  n  will  be  different  from  the  value  of  the  electric  field  perpendicular  to  the 


current  Eo1  .  The  probability  p(&)  of  finding  the  resulting  electric  field  vector  E0  in  an 
interval  d0  around  0  will  then  correspond  to  an  ellipsoid  given  by 


p(Q)d6 


Vcos2  0  +z4  sin2  0  sin QdQ 
(cos2  0  +  z2  sin2  0) 


(2.45) 


where 

z  =  EoU  /EoL=a!b,  (2.46) 

where  a  is  along  the  direction  of  the  current,  and  b  is  perpendicular  to  the  direction  of  the 
current,  as  illustrated  in  Figure  2-3.  We  can  then  obtain  the  profile  for  the  ion  being 
subject  to  a  large  number  of  waves  with  different  phase  by  applying  the  Rayleigh 
distribution  to  Eq.  (2.38) 


lm  =<  Inm  >  £„  *  |(' V  „  K*)|  V«f  f  0  M))  £o  5[A®  +  P®]  .  (2-47) 

p—~c o 


where 


(j2r(^E0Ak))  =jy’(4£,M) 
0  0 


(2.48) 
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Letting  x  =  QE0  Ak  we  can  write 


(2.49) 


(2.50) 


where 


Q iEpAkf  (sM)! 
2  2 


(2.51) 


The  parameter  s  in  the  model  corresponds  to  the  value  of  the  electric  field  strength  of  the 
waves  and  is  given  by 


3  hE0 

e  = - - — . 

2eme(opeZ 


(2.52) 


As  discussed  in  Section  2.1,  the  source  of  oscillating  electric  fields  in  the  Phoenix  plasma 
is  assumed  to  be  Langmuir  waves,  so  co  has  been  replaced  by  the  plasma  frequency 


l4%Ne. 

i  me 


(2.53) 


The  line  profile  for  a  single  plasma  frequency  corresponding  to  the  plasma  wavelength 

l2'54) 

will  then  be  given,  apart  from  the  dipole  term,  by (1) 
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(2.55) 


For  Efc  »  1 ,  the  envelope  of  the  satellites  is  Gaussian  (Figure  2-4)  with  a  width 


A.  [A„ (cm)]2  a0e3 eCV2 Ak  , - T. - -r-  - 

^ hwhm (cm)  ~  i  ~2  ^471(1.56  *10  cm )m  , 

2ncn 


(2.56) 


which  when  simplified  and  put  in  terms  of  mA  becomes 


a?WaA™ A)  *  3.72  •  10_7(?b0(mA))2  sC1/2M . 


(2.57) 


For  Lp  with  eC  Ak  =  1,  this  results  in  a  HWHM  of  13.6  mA. 


2.2.3.2  Power  Spectrum  of  cope  Due  to  collisional  damping  and  variations  in  density, 
there  will  be  a  distribution  of  oscillating  electric  field  frequencies.  The  model  includes 
this  by  using  a  Lorentz  distribution  for  the  plasma  frequencies,  that  is, 


L(Aco,g)  = 


2k  a  2  (g^pe 

Aar  +  — -- 

l  2  . 


(2.58) 


where  the  corresponding  HWHM  of  this  distribution  is 


g  =  g^Pe  g  47ie2  (1.56 -1022  cm'3  )C 
2  2  V  m„ 


(2.59) 


6(rad /  s)  =  3.5  x  10l5gC1/2 . 


(2.60) 
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Intensity 


Waveiength  (mA) 


Figure  2-4.  Splitting  of  a  line  due  to  an  oscillating  electric  field. 
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The  parameter  g  characterizes  the  width  of  the  power  spectrum  of  electric  field 
oscillations. 

The  line  profile  due  to  the  oscillating  electric  fields  with  this  distribution  of  applied 
frequencies  is 

'(A»-)=  I /,(s)i(AA-,A>!s).  (2.61) 

p=- CO 

The  width  of  the  Lorenz  profiles  depend  upon  the  satellite  number  \p\ .  As  \p\  increases, 
the  width  of  the  Lorentz  profiles  also  increase. 

2.2.3  SELF  ABSORPTION. 

Radiation  emitted  by  ions  must  make  it  out  of  the  plasma  before  it  can  be  detected  by  a 
spectrograph.  Photons  that  are  emitted  by  the  ion  may  be  absorbed  by  other  ions  before 
they  escape  the  plasma.  The  newly  excited  atoms  from  these  photons  can  re-emit  the 
photon  or  lose  the  energy  by  collisional  de-excitation  resulting  in  less  photons  escaping 
the  plasma  than  were  emitted.  Radiation  from  the  line  center  will  most  strongly  be 
absorbed,  since  at  line  center,  where  there  are  many  atoms  emitting,  there  are  a  large 
number  of  atoms  in  the  correct  state  to  absorb  a  photon.  As  a  result,  self  absorption 
reduces  the  line  center  more  than  the  wings.  In  the  model  we  use  the  standard  expression 
for  self  absorption/23 ) 
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f 


(2.62) 


/(AX.)  = 


l—e  /ro 


where 

t  =  the  optical  depth  of  the  spectral  line  at  the  maximum  of  the  profile  /a(AX.). 


2.2.4  ZEEMAN  BROADENING. 


We  must  also  consider  the  possibility  that  Zeeman  splitting  is  the  principal  broadening 
mechanism  in  the  plasma/10'1 2)  The  difference  in  the  Zeeman  splitting  of  the  levels  could 
explain  the  line  dependent  differences  in  the  observed  profiles.  We  consider  this 
possibility  by  determining  the  magnetic  fields  required  to  produce  broadening  on  the 
order  of  the  observed  lines. 

For  strong  fields,  the  change  in  frequency  associated  with  the  Zeeman  splitting  of  an 
atomic  level  will  be  given  by 


A  ca  = 


VbB(T) 


(2.63) 


which  corresponds  to  a  change  in  wavelength  of (14) 

A X(m)  =  Aco  =  Xl^bB{T)  .  (2.64) 

2nc  Inch 

Inserting  the  appropriate  constants  and  converting  to  Gaussian  units  we  arrive  at  an 
expression  for  the  shift  in  wavelength  of  a  line  for  a  given  wavelength  X0  and  magnetic 
field  B,  namely. 
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AZ(mA)  =  4.63  x  1  O'10  \0(A)]2  B(G ) . 


(2.65) 


To  produce  the  width  of  the  observed  profiles  for  the  Lyman  series,  which  are  on  the 
order  of  30  mA,  a  magnetic  field  of  2.0  gigagauss  (GG)  would  be  required.  This  would 
result  in  a  ratio  of  the  magnetic  to  kinetic  energy  density  of  100  for  a  plasma  temperature 
of  T  «  1  keV.  From  the  best  available  models  of  z-pinches,  this  ratio  will  rarely  exceed 
unity.  Physically,  the  current  tends  to  filament  around  the  pinch  resulting  in  low 
magnetic  fields  inside  the  pinch.  Zeeman  splitting  cannot  be  the  principle  broadening 
mechanism  for  the  profiles. 


2.3.  COMPLETE  LINESHAPE  FOR  THE  PHOENIX  PLASMA  SOURCE. 

The  spectral  profile  for  the  Phoenix  plasma  source  without  self  absorption  is  calculated 
by  convoluting  the  Doppler  profile  Eq.  (2.5),  the  Holtzmark  profile  Eq.  (2.14),  and  the 
oscillating  electric  field  profile  Eq.  (2.61).  By  using  the  results  of  the  convolution  of  a 
Doppler,  Holtzmark,  and  Lorentzian  profile,  as  derived  in  Appendix  A, 


D(AZ,a,p,5)  =  —  \ye 

rr  J 
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[AZ  sin(PAAy) + §  cos(PAXy)  +  cos(pAXy)}fy ,  (2.66) 
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the  profile  for  each  Stark  component  without  self  absorption  will  be  given  by 


/(AX)=  EV  zD[AX-p\p,a$,\p\g).  (2.67) 

p—— X 
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The  final  profile  is  calculated  by  summing  up  the  profiles  for  each  Stark  component  and 
applying  the  expression  for  self  absorption. 


2.4.  SUMMARY  OF  PLASMA  PARAMETERS. 

We  have  presented  here  a  model  which  attempts  to  take  into  account  all  of  the  broadening 
mechanisms  present  in  dense  plasmas,  such  as  produced  on  the  Phoenix  device.  The 
parameters  introduced  for  the  model  are:  b,  C,  z,  s,  g  and  t.  Each  parameter  gives  some 
information  about  the  plasma. 

The  parameter  b  is  a  measure  of  the  velocity  of  the  bulk  motion  of  the  plasma.  It  will 
make  its  main  contribution  in  the  center  of  the  line.  For  the  central  component  of  Ly  and 
Ls  it  will  be  the  only  broadening  mechanism. 

The  parameter  C  is  a  measure  of  the  density  of  the  plasma.  This  parameter  effects  the 
line  shape  through  all  forms  of  Stark  broadening.  The  Stark  broadening  by  ions  will 
make  large  contributions  to  the  wings  of  the  line  since  the  Holtzmark  profiles  fall  off 
slower  than  Gauss  and  Lorentz  profiles.  C  also  determines  the  frequency  of  the 
oscillating  electric  fields  through  the  plasma  frequency  and  C  makes  a  small  contribution 
the  width  of  the  envelope  of  the  oscillating  electric  field  profile. 

The  parameter  z  is  a  measure  of  the  degree  of  anisotropy  of  the  oscillating  electric  fields. 
A  value  of  z  =  1  corresponds  an  equal  probability  of  finding  the  direction  of  propagation 
in  all  directions.  For  z  >  1 ,  the  probability  of  finding  the  propagation  along  the  pinch  axis 
is  higher  than  perpendicular  to  the  pinch  axis. 
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The  parameter  s  measures  the  magnitude  of  the  oscillating  electric  fields.  This  field 
strength  variable  determines  the  width  of  envelope  of  the  satellite  lines  under  the  action 
oscillating  electric  fields. 

The  parameter  g  is  a  measure  of  the  width  of  the  power  spectrum  of  the  oscillations  in  the 
plasma.  Since  the  frequency  of  the  oscillations  depends  on  density,  this  value  also  gives 
some  measure  of  density  fluctuations. 

And  finally,  the  parameter  t  is  the  optical  depth  at  the  center  of  the  line  profiles.  This  is  a 
measure  of  the  absorption  of  the  radiation  before  it  makes  it  out  of  the  plasma. 


SECTION  3 

EXPERIMENTAL  APPARATUS 


3.1  PHOENIX  DEVICE. 

Phoenix  PRS  is  a  pulsed  power  z-pinch  device  used  for  producing  radiation  in  the  soft  x- 
ray  region.  A  diagram  of  the  machine  is  shown  in  Figure  3-1 .  The  two  primary 
components  of  the  machine  are  the  pulser  and  the  load.  The  pulser  generates  a  fast  rise 
time  power  pulse  to  deliver  to  the  load,  which  typically  consists  of  a  wire  array  or  a  gas 
puffed  into  the  apparatus. 

The  schematic  circuit  diagrams  of  the  Phoenix  machine  are  shown  in  Figure  3-2.(26)  The 
analysis  of  the  power  flow  in  this  device  is  quite  complex.  To  simplify  the  problem,  the 
machine  is  often  modeled  in  terms  of  an  equivalent  circuit  as  in  Figure  3-3.  By  knowing 
how  the  machine  couples  to  the  load,  the  reflected  power  from  the  load  can  be  minimized 
for  optimal  performance.  For  the  data  in  this  dissertation,  the  load  consisted  of  sixteen  2 
cm  aluminum  wires,  0.8  mil  in  diameter,  arranged  in  a  circular  array  of  diameter  2.6  cm 
(Figure  3-4).  The  machine  typically  delivers  a  current  of  4  Ma  to  the  load  and  produces 
10-50  kJ  of  energy  in  the  Lyman  radiation  of  aluminum  in  this  configuration. 

3.2  DIAGNOSTICS. 

The  diagnostics  used  to  study  the  plasma  produced  by  Phoenix  were  a  PIN  diode 
detector,  a  diamond  PCD,  a  flat  crystal  spectrometer,  a  pinhole  camera  and  two  bent 
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Figure  3-1.  Phoenix  advanced  radiation  source. 
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Figure  3-2.  Schematic  diagram  of  the  Phoenix  pulsed  radiation 
source. 
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Firuge  3-3.  Equivalent  circuit  for  the  Phoenix 
pulsed  radiation  source. 
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crystal  spectrographs  of  the  Johann  type.  The  location  of  the  diagnostics  on  the  Phoenix 
spool  are  shown  in  Figure  3-5. 

3.2.1  X  Ray  Detectors. 

The  radiation  pulse  from  the  Phoenix  device  was  monitored  using  a  diamond 
photoconductive  detector  (PCD)  and  a  Quantrad  PIN  diode.  The  detectors  give  the 
intensity  of  the  radiation  as  a  function  of  time  for  the  radiation  event.  The  total  energy 
can  be  obtained  from  the  integration  of  the  signals  obtained  from  these  devices.  The 
diamond  PCD  consists  of  a  thin  diamond  crystal  and  the  PIN  diode  consists  of  an  ultra 
high  purity  silicon  (nearly  Intrinsic)  wafer,  with  one  side  doped  with  a  P-layer  and  the 
other  side  doped  with  an  N-layer  leading  to  the  so  called  PIN  diode  structure.  Both 
detectors  are  operated  by  applying  a  bias  voltage  across  the  active  material.  When  an  x 
ray  undergoes  a  collision  with  an  electron,  an  electron-hole  pair  is  created,  and  the 
electrons  resulting  from  these  collisions  are  swept  out  by  the  bias  voltage.  The  PIN  diode 
is  operated  with  a  negative  bias  voltage  and  has  a  sensitivity  about  1000  times  greater 
than  the  diamond  PCD. 

The  response  curves  for  these  devices  are  shown  in  Figure  3-6.  The  response  curve  of  the 
diamond  PCD  was  provided  by  the  manufacturer  at  Los  Alamos  National  Laboratory. 

The  response  curve  for  the  PIN  diode  was  calculated  using  a  formula  supplied  by  the 
manufacturer  (Quantrad). 

Since  we  are  interested  in  the  most  energetic  x  rays,  corresponding  to  the  highest 
ionization  state  of  the  atoms,  filters  were  placed  in  front  of  the  detectors  to  strongly 
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attenuate  the  longer  wavelengths.  The  filters  also  served  to  prevent  debris  from 


damaging  the  detectors  and  to  prevent  the  detectors  from  becoming  saturated.  The  filters 
used  with  the  diamond  PCD  detector  were  0.3  mil  Kapton,  0.16  mil  Kimfoil,  and  3000  A 
aluminum.  A  400  p  m  diameter  pinhole  in  a  2-mil  stainless  steel  substrate  was  used  with 
the  PIN  diode  to  prevent  the  detector  from  becoming  saturated.  The  pinhole  was 
mounted  in  a  2  mm  diameter  hole  placed  on  a  8  mm  thick  brass  disk  and  covered  with  a 
5-mil  Be  foil  as  shown  in  Figure  3-7. 

The  absorption  curves  for  the  filters  used  in  the  detectors  are  shown  in  Figure  3-8.  The 
mass  absorption  coefficients  for  the  filters  were  computed  using  the  empirical  formula 
given  by  Victoreen 

^C),3-DX4  +  Pk-nN?-  (3.1) 

P  A 


with  X  the  wavelength,  aK_N  the  Klein-Nishina  cross  section,  TV  A vogadros  number,  Z 
the  atomic  number,  and  A  the  atomic  Weight.'275  The  Klein-Nishina  cross  section  is 
given  by 


N  —  2  7Z 


1  4 

me 


1  +  a  [2(1 -i- a)  1  ,  _  .1  1  W1  _  .  1  +  3a 

— - ln(l+2a)[  +  —  ln(l+2a)-- — — r 

or  1  2a +  1  a  j  2a  (l+2a)_ 


(3.2) 


where  a  =  hv/  me 2  .  Since  the  actual  measured  values  of  the  mass  absorption  coefficients 
are  intermittent  in  the  wavelength  regions  of  interest,  between  5  A  and  20  A,  it  was 
necessary  to  use  the  empirical  formula  given  by  Eq.  (3.1).  The  measured  absorption 
coefficients  were  fit  to  find  C  and  D.(28;  Since  the  absorption  coefficients  change 
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(a)  PIN  diode  detector. 
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(b)  Pinhole  camera. 


Figure  3-7.  Arrangement  for  PIN  diode  and  pinhole  camera 
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Figure  3*8.  X  ray  absorption  curves. 
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abruptly  at  the  absorption  edges  of  the  atom,  different  values  for  C  and  D  were  computed 
for  the  different  regions  of  the  spectrum.  The  total  response  of  the  detectors  was 
computed  from  the  product  of  the  response  of  the  detectors  with  their  filters  as  shown  in 
Figure  3-9.  The  respective  distances  of  the  Diamond  PCD  and  the  PIN  Diode  from  the 
wire  array  were  2.92  m,  and  1 .69  m. 

3.2.2  Pinhole  Camera. 

A  pinhole  camera  as  shown  in  Figure  3-7  was  utilized  to  get  time  integrated  images  of  the 
x  rays  emitted  from  the  plasma.  The  pinhole  pictures  measure  the  spatial  distribution  of 
the  x-ray  intensity  in  the  source.  Two  configurations  of  the  camera  were  used,  but  for 
both  configurations  the  distance  from  the  source  to  the  pinhole  was  169  cm.  For  shots 
1118-1 125,  the  film  was  placed  50  cm  from  a  400  pm  pinhole  resulting  in  a 
magnification  of  0.3.  The  filters  used  for  this  set  of  shots  consisted  of  2  mil  A1  with  5  mil 
Be.  For  shots  1 126-1 129,  the  film  was  placed  69  cm  from  a  25  pm  pinhole  resulting  in  a 
magnification  of  0.4.  The  filters  for  these  shots  consisted  only  of  5  mil  Be. 

KodakBiomax  MR  film  was  used  because  of  its  high  sensitivity  and  clear  background 
emulsion.  The  absorption  curves  for  the  pinhole  filters  are  given  in  Figure  3-8. 

3.2.3  Four  Frame  Camera. 

A  four  frame  camera  was  used  to  determine  the  size  of  the  pinch  events  and  to  observe 
the  temporal  evolution  of  the  plasma.  Time  resolved  imaging  is  necessary  to  determine 
the  minimum  plasma  column  size,  since  in  time  integrated  imaging  the  final  size  of  the 
pinch  is  masked  by  images  from  other  parts  of  the  plasma  at  other  times.  The  major 
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Wavelength  (A) 

(a)  Diamond  PCD  response  with  2  mil  Kapton  and  3000  A  Al. 


Wavelength  (A) 

(b)  PIN  diode  response  with  5  mil  Be. 


Figure  3-9.  Comparison  of  PCD  detector  response  curve  versus  PIN  diode. 


components  of  the  four  frame  camera  are  an  MCP,  a  phosphor  screen,  and  a  fiber  optic 
array  as  illustrated  in  Figure  3-10.(29)  The  MCP  consists  of  10-100  pm  diameter  glass 
capillary  tubes  stacked  side  to  side.  When  a  voltage  is  applied  across  the  ends  of  the 
tubes,  the  tubes  perform  as  do  photomultiplier  tubes,  shown  in  Figure  3-10.  The  resulting 
cascade  of  electrons  impact  a  phosphor  screen  where  they  cause  the  phosphor  to  emit 
photons  in  the  visible  region  of  the  spectrum.  Since  the  lifetime  of  the  phosphor  is  long, 
typically  40-80  ns,  an  image  can  be  obtained  on  standard  film. 

The  “on/off  ’  state  of  the  MCP  is  controlled  by  the  voltage  applied  across  it.  If  a  voltage 
pulse  is  applied  to  the  MCP,  the  camera  is  ON  only  for  the  duration  of  this  pulse,  since 
the  gain  of  the  device  is  negligible  when  no  voltage  is  applied,  and  it  is  approximately  106 
when  the  optimal  voltage  is  applied.  In  the  4FC  the  MPC  is  cut  into  four  sections  which 
are  electrically  isolated  from  one  another  allowing  the  four  sections  of  the  camera  to 
operated  individually  as  shown  in  Figure  3-1 1.  Applying  short  voltage  pulses,  typically  1 
to  10  ns,  to  each  section  of  the  camera  sets  the  exposure  time  for  each  frame.  The  pulses 
sent  to  each  section  of  the  camera  are  delayed  by  1  to  30  ns  between  frames  resulting  in 
framing  of  the  camera.  For  the  experiments  performed  at  Phoenix,  the  frames  were 
opened  for  5  ns  with  a  5  ns  delay  between  frames. 

Figure  3-1 1  shows  the  setup  used  to  image  a  plasma  with  the  4FC.  Pinholes 
are  placed  in  front  of  each  section  of  the  4FC  to  produce  images  on  the  corresponding 
sections  of  the  MCP.  For  the  experiment,  the  distance  from  the  pinhole  to  the  source  was 
1 70  cm  and  the  distance  from  the  pinhole  to  the  film  was  22  cm,  resulting  in  a 
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magnification  of  0.13.  Pinholes  of  100  pm  and  400  pm  were  used  to  ensure  that  enough 
intensity  was  reaching  the  camera  in  the  short  time  that  the  cameras  were  on.  With  such 
large  pinholes  compared  to  the  wavelength,  the  resolution  of  the  camera  was  limited  by 
the  size  of  the  pinholes. 

An  important  aspect  to  the  operation  of  the  4FC  is  the  timing  of  the  camera.  Ideally  the 
camera  should  start  framing  just  before  the  pinch  event  occurs  so  that  the  formation  of  the 
pinch  can  be  observed.  Figure  3-12  shows  the  timing  scheme  used  to  open  the  camera  at 
the  correct  time.  The  timing  begins  when  a  signal  from  a  V-dot  probe,  a  voltage  monitor 
which  detects  the  pulse  in  the  machine  250  ns  prior  to  the  pinch,  triggers  the  delay 
generator.  The  delay  generator  delays  the  pulse  about  200  ns  and  triggers  the  four  frame 
camera  pulser.  The  pulser  sends  high  voltage  pulses  to  the  four  frames  of  the  camera 
with  respective  delays  to  each  frame.  The  pulse  returning  from  the  first  frame  is  then  sent 
to  the  control  room.  The  radiation  from  the  machine  is  monitored  using  the  PIN  diode, 
described  in  Section  3.2.1,  and  this  signal  is  correlated  to  the  pulse  which  returned  from 
the  4FC.  Adjustments  of  relative  timing  are  made  using  a  Stanford  Research  Systems 
model  DG535  delay  generator  to  allow  for  proper  placement  of  the  4FC  triggering  in  the 
radiation  pulse. 

3.2.4  Spectrographs. 

3.2.4.1  Flat  Crystal  Spectrograph.  The  flat  crystal  spectrograph  was  employed  to  look 
at  lines  in  the  Lyman  series  of  A1 XIII.  In  particular  the  spectrograph  was  designed  to 
detect  the  Ly,  Lg  and  Ls  lines.  The  data  obtained  from  this  spectrometer  was  used  to 
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Figure  3-12.  Four  Frame  Camera  Timing. 
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check  the  results  of  the  bent  crystal  spectrographs  and  give  additional  information  about 
parameters  in  the  plasma.  The  Auburn  flat  crystal  spectrometer  shown  in  Figure  3-1 3(a) 
uses  two  flat  crystals  instead  of  one,  since  it  is  a  modification  of  a  previous  diagnostic. 
Two  Quartz  1 0 J_0  crystals  of  3. 12  cm  were  used  in  the  Auburn  spectrograph.  The  crystals 
were  set  for  an  angle  of  a  «  42.5°.  The  source  to  crystal  distance  was  U  =  100  cm  and 
the  crystal  to  film  distance  was  V  =  10  cm.  This  configuration  results  in  a  detectable 
wavelength  region  of  5.61  A  to  5.89  A  which  encompasses  the  Ly ,  L~ ,  and  Lt  lines 
having  central  wavelengths  of  5.739  A,  5.605  A,  and  5.534  A  respectively.  From 
geometrical  considerations,  the  wavelength  as  a  function  of  position  was  calculated  for 
the  spectrograph  (Appendix  B).  The  dispersion  is  linear  for  these  spectrograph 
parameters  and  the  value  of  the  dispersion  calculated  was  5.8  mA/mm.  The  angle  a 
could  be  set  for  each  shot  to  within  0.5°,  however,  vibrations  in  the  machine  contributed 
to  fluctuations  in  a  from  shot  to  shot.  A  more  accurate  value  of  a  was  later  calculated 
from  the  spectra  using  the  relative  positions  of  two  known  lines.  A  5  cm  long,  1  mm 
wide  slit  was  placed  between  the  source  and  crystal  located  12  cm  from  the  crystal.  The 
slit  was  oriented  along  the  direction  of  dispersion  of  the  crystal  as  shown  in  Figure  3-13. 
A  2-mil  Be  foil  was  placed  over  the  slit  to  filter  out  longer  wavelengths  and  prevent 
saturation  of  the  film. 

The  resolution  of  flat  crystal  spectrometers  is  generally  limited  by  the  quality  of  the 
crystal  used  and  the  spatial  extension  of  the  source.^30'3  ^  The  quality  of  the  crystal  affects 
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(b)  Instrumental  broadening  due  to  geometry. 


Figure  3-13.  Schematic  of  flat  crystal  spectrometer. 
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the  resolution  through  the  variation  of  the  2d  spacing  in  the  crystal  and  through  the 
surface  structure.  The  2d  value  for  real  crystals  will  vary  slightly  due  to  imperfections  in 
the  crystals.  This  variation  of  the  2d  value  broadens  the  spectra  on  the  detector.  Using  the 
resolving  power  for  quartz  1010  of  X  /  A X  =  7.54  x  103,  the  resulting  resolution  is  2  mA 
for  the  Auburn  flat  crystal  spectrograph.  The  effect  of  the  spatial  distribution  of  the 
source  on  the  resolution  arises  from  the  geometry  of  the  spectrograph.  For  a  flat  crystal 
spectrograph  this  broadening  can  be  quite  large  depending  upon  the  size  of  the  source. 
This  is  illustrated  in  Figure  3-13.  We  consider  two  points  in  the  source  separated  by  a 
distance  h  in  the  direction  of  dispersion  of  the  crystal,  and  we  look  at  two  parallel  rays 
from  different  parts  of  the  source  that  have  the  same  wavelength  and  satisfy  the  Bragg 
angle  a.  The  rays  from  these  two  points  with  the  same  wavelengths  will  strike  different 
places  on  the  film.  This  results  in  the  broadening  of  the  line.  There  are  a  number  of  ways 
for  two  rays  of  the  same  wavelength  from  different  parts  of  the  plasma  to  strike  the 
crystal  at  the  same  Bragg  angle.  The  sum  of  these  effects  can  be  calculated  using  ray 
tracing.  The  geometry  of  the  Auburn  flat  crystal  spectrometer  gives  a  displacement  of 
1  mm  on  the  film  when  there  is  a  displacement  of  1  mm  in  the  source.  From  the 
previously  calculated  dispersion  of  5.8  mA/mm,  this  means  that  for  every  mm  that  the 
source  is  extended  along  the  direction  of  dispersion  of  the  crystals,  the  line  is  broadened 
by  5.8  mA.  For  this  reason  the  spectrograph  was  oriented  so  that  the  smaller  dimension 
of  the  plasma  was  along  the  dispersion  direction  of  the  crystals.  As  we  will  see  in  the 
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next  section,  curved  crystal  spectrometers  in  the  Johann  configuration  do  not  always  have 
a  large  broadening  due  to  the  geometry. 

3.2.4.2  Bent  Crystal  Spectrographs.  The  primaiy  diagnostics  used  for  measuring  the 
polarization  lineshapes  in  the  plasma  were  two  bent  crystal  spectrometers  mounted  in  the 
Johann  configuration  (Figure  D-l).  These  two  spectrometers  were  used  to  obtain  high 
resolution  line  profiles  of  the  and  L.(  lines  of  A1 XIII  in  two  states  of  polarization. 

On  each  shot  both  spectrometers  monitored  the  same  line,  either  Zp  or  Ly . 

Both  Johann  spectrographs  used  a  quartz  1010  crystal  bent  to  a  radius  of  curvature  of 

1- m.  A  5  cm  long,  1  mm  wide  slit  was  placed  between  the  source  and  crystal,  1 1  cm 
from  the  crystal.  The  slit  was  oriented  parallel  to  the  plane  of  the  spectrograph  and  a 

2- mil  Be  filter  (Figure  3-8)  was  placed  over  the  slit. 

The  crystals  were  oriented  so  that  the  plane  of  polarization  of  one  of  the 
spectrometers  was  parallel  to  the  bulk  current  in  the  plasma,  and  the  plane  of  polarization 
of  the  other  was  perpendicular  to  the  bulk  current  in  the  plasma  as  illustrated  in  Figure  3- 
14.  The  details  of  how  x  rays  are  polarized  by  crystals  under  Bragg  reflection  are  given 
in  Appendix  C.  The  difference  in  these  two  configurations  results  in  one  crystal  being 
illuminated  by  a  point  source  along  the  dispersion  of  the  crystal  and  the  other  being 
illuminated  by  a  line  source.  To  see  how  the  orientation  of  the  source  would  affect  the 
observed  profiles,  we  performed  ray  tracing  for  the  Johann  spectrometers.  The  details  of 
the  ray  tracing  calculations  for  the  Johann  spectrometers  is  given  in  Appendix  D. 
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Crystal  A 


Crystal  B 


3.2.4.2.1  Ray  Tracing  Results  for  the  Bent  Crystal  Spectrographs.  For  the  ray 
tracing  calculations,  the  source  to  crystal  distances  tested  were  1.5  m  and  3.0  m.  These 
two  distances  were  chosen  based  upon  required  intensities  at  the  film,  and  this  was  done 
for  two  crystal  lengths  of  2.5  cm  and  7.5  cm.  The  input  sources  for  the  calculations  were 
a  point  source  and  a  2  cm  long  line  source.  Ray  tracing  was  performed  for  wavelengths 
of  interest  for  the  Phoenix  experiments  of  6.053  A  and  5.739  A  corresponding  to  the  Z,p 
and  Ly  lines  of  A1 XIII.  Gaussian  profiles  with  constant  radiance  centered  at  these 

wavelengths  were  input  with  different  half  widths  for  the  sources  to  test  the  limitations  of 
the  spectrographs  in  various  configurations.  For  a  well  designed  spectrograph,  the  output 
spectra  should  reconstruct  the  profile  in  the  source. 

Before  ray  tracing  was  performed  using  Gaussian  profiles  in  the  source,  the  resolution 
and  detectable  wavelength  of  the  test  configurations  were  found.  The  resolution 
determines  the  minimum  difference  in  wavelength  that  can  be  resolved  by  the 
spectrometer  and  the  detectable  wavelength  region  determines  the  maximum  width  of 
lines  that  can  be  observed. 

The  resolution  of  the  spectrographs  for  the  sources  are  determined  by  ray  tracing  a  single 
wavelength  for  every  point  in  the  source/30  By  putting  in  an  infinitely  narrow  line  at  a 
single  wavelength  and  observing  the  width  of  the  output  line  we  determined  the 
broadening  due  to  the  instrument.  The  results  of  these  calculations  for  the  and  Ly 
lines  of  A1  XIII  are  shown  in  Figure  3-15.  The  resolving  power  for  the  point  source, 
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AX,  (mA) 

(a)  Resolving  power  for  a  point  source 


AX,  (mA) 

(b)  2  cm  line  source 


Figure  3-15.  Intensity  versus  wavelength  for  point  and  extended  sources 


shown  in  Figure  3-15,  obtained  from  this  ray  tracing  technique  is  approximately  0.82 
mA  for  Ip  and  0.77  mA  for  Ly ,  independent  of  the  source  to  crystal  distance.  For  the 
line  source,  the  resolution  depends  on  both  the  wavelength  and  distance  of  the  crystal 
from  the  source  as  shown  in  Figure  3-15.  For  Ip ,  the  resolving  power  for  the  line  source 
is  approximately  0.89  mA  at  a  distance  of  1 .5  m  and  0.82  mA  at  3.0  m.  For  L ,  the 

resolving  power  for  the  line  source  is  approximately  0.86  mA  at  1.5  m  and  0.78  mA  at 
3.0  m.  In  the  Johann  configuration,  the  primary  factors  effecting  the  resolution  are  the 
diameter  of  the  Rowland  circle  and  the  quality  of  the  crystal. 

The  detectable  wavelength  region  of  the  spectrographs  were  found  by  calculating  the 
angle  that  the  rays  from  the  source  points  make  with  the  ends  of  the  crystal  and  then 
converting  these  angles  to  wavelength  through  Bragg’s  law.  The  detectable  wavelength 
regions  for  the  Lp  and  Ly  lines  for  three  points  on  a  line  source  are  plotted  for  crystal 
lengths  of  2.5  cm  in  Figure  3-16  and  7.5  cm  in  Figure  3-17.  The  difference  in  the 
detectable  wavelength  region  due  to  the  geometry  is  illustrated  by  these  plots.  For  the  L? 
line  and  Ly  line,  the  detectable  wavelengths  regions  are  approximately  the  same  because 

the  geometry  of  the  spectrograph  does  not  change  much  between  these  lines.  In  designing 
a  spectrograph,  the  detectable  wavelength  should  be  greater  than  twice  as  large  as  the  half 
width  of  the  observed  line  in  order  to  observe  the  wings  of  the  profiles.  Since  the  largest 
lines  expected  for  the  Phoenix  experiment  are  on  the  order  of  50  mA  FWHM,  a  crystal 


55 


Wavelength  (A) 


Distance  (m) 
(a)  A  =  5.739  A 


Distance  (m) 
(b)  A  =  6.053  A 


Figure  3-16.  Detectable  wavelength  region  for  a  Johann  spectrometer  with  a 
2.5  cm  crystal. 
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Distance  (m) 
(b)  k  =  6.053  A 


Figure  3-17.  Detectable  wavelength  region  for  a  Johann  spectrometer  with  a 
7.5  cm  crystal. 
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length  of  2.5  cm  is  not  sufficient.  The  plots  of  detectable  wavelength  show  that  the  7.5 
cm  crystal  should  be  used  for  the  experiments  on  Phoenix. 

Test  Gaussian  spectral  profiles  were  then  input  for  the  point  and  line  sources  and 
the  output  profiles  of  the  spectrographs  were  compared  to  the  input  profiles.  Gaussian 
profiles  of  from  20  mA  to  1 50  mA  FWHM  were  traced  for  the  test  configurations.  The 
output  lines  for  a  given  input  half  width  are  shown  in  Figures  3-18  through  3-21  for  the 
2.5  cm  and  7.5  cm  crystals.  The  results  for  the  L.,  line  are  not  presented  because  they  do 
not  vary  much  from  the  results  for  the  ip  line.  These  plots  show  that  for  a  line  of  half 
width  of  50  mA,  the  profile  for  both  a  point  and  2  cm  line  will  be  correctly  measured  by 
the  spectrograph  if  the  7.5  cm  crystal  is  used. 

The  data  gained  from  these  calculations  was  used  to  design  a  spectrometer  which  would 
perform  properly  for  Phoenix  plasma  source.  The  three  factors  of  primary  concern  for 
designing  the  spectrograph  were  the  source  to  crystal  distance,  the  crystal  length  and 
resolution.  The  distance  of  the  source  to  the  spectrometer  was  important  for  achieving 
the  correct  intensity  at  the  film.  From  the  graphs  of  Figure  3-16  and  3-17  we  see  that  the 
detectable  wavelength  region  is  proportional  to  the  size  of  the  crystal  and  the  source  to 
crystal  distance.  The  usual  procedure  is  to  choose  the  diameter  of  the  Rowland  circle 
required  to  get  the  desired  resolution,  find  the  distance  which  gives  the  correct  intensity 
and  then  use  a  crystal  of  the  appropriate  size.  If  this  approach  fails  then  the  diameter  of 
the  Rowland  circle  must  be  changed  to  resolve  the  problem.  From  the  ray  tracing  it  was 
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Wavelength  (A) 

(a)  7.5  cm  for  a  point  source 


Wavelength  (A) 
(b)  2  cm  line  source 


Figure  3-19.  Ray  Tracing  for  Lp  with  a  source  to  crystal  distance  of  1.5  m  and 
a  crystal  length. 
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a) 


b) 


Figure  3-20.  Ray  Tracing  for  Lp  with  a  source  to  crystal  distance  of  3.0  m  and 
a  crystal  length  of  2.5. 
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determined  that  the  best  configuration  for  the  Phoenix  machine  was  a  crystal  length  of 
7.5  cm.  and  a  source  distance  of  1.5  m. 


63 


SECTION  4 

EXPERIMENTAL  RESULTS 


4.1  RADIATION  DETECTOR  RESULTS. 

The  radiation  detectors  provided  information  about  the  time  evolution  of  the  radiation 
from  the  plasma.  Typical  signals  obtained  from  the  Si  PIN  diode  and  diamond  PCD  are 
shown  in  Figure  4-1 .  Temporally  the  detectors  are  in  agreement  about  the  time  of  arrival 
of  the  radiation  pulse,  however,  the  intensity  seen  on  the  two  detectors  shows  a 
considerable  difference  in  the  tail  of  the  pulse.  The  PIN  diode  detector  has  a  larger  signal 
after  the  primary  pulse.  The  difference  in  the  signals  from  the  two  detectors  cannot  be 
explained  by  the  filters  and  response  curves  of  the  detectors  alone  (Section  2.1). 
According  to  the  calculations  for  the  response  of  the  detectors,  the  Si  PIN  diode  will  have 
a  lower  wavelength  cutoff  than  the  diamond  PCD.  This  would  suggest  that  the  PIN  diode 
should  have  a  shorter  pulse,  since  the  lower  wavelengths  should  occur  for  a  smaller 
duration.  This  difference  in  the  two  signals  can  only  be  accounted  for  by  the  sensitivity 
of  the  Si  PIN  diode  to  visible  light,  for  the  epoxy  used  to  seal  the  Be  foils  was  found  to  be 
partially  transparent  for  high  intensities  of  light. 

The  other  useful  information  provided  by  these  detectors  is  the  power  output  in  the 
radiation.  The  total  energy  of  the  radiation  pulse  for  each  shot  was  calculated  by 
integrating  the  signal  from  the  PCD  detector.  The  total  energy  in  the  radiation,  for  all 
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wavelengths  below  10  A,  from  the  PCD  is  given  for  each  shot  in  Table  4-1,  along  with 
the  peak  machine  current. 

Table  4-1 .  Peak  current  and  total  energy  output  of  the  phoenix  device. 


Phoenix  Shot  Number 

Peak  Current  (MA) 

Total  Output  Energy  (kJ) 

1109 

3.3 

14.5 

1110 

2.9 

10.5 

1111 

3.5 

19.4 

1112 

3.4 

16.6 

1113 

3.3 

16.0 

1114 

3.5 

20.9 

1115 

3.4 

15.7 

1116 

3.5 

19.7 

1117 

3.6 

16.7 

1118 

3.5 

23.3 

1119 

3.4 

17.2 

1120 

3.4 

19.0 

1121 

3.4 

15.8 

1122 

3.5 

15.2 

1123 

3.4 

17.6 

1124 

3.4 

14.5 

1125 

3.5 

19.9 

1126 

3.5 

20.2 

1127 

3.5 

21.2 

1128 

3.5 

21.8 

1129 

3.5 

19.6 

Table  4-1  shows  that  there  is  not  much  variation  of  the  peak  current  even  though  the 
radiation  output  varies  from  shot  to  shot.  It  is  conceivable  that  the  variation  in  the  output 
radiation  results  from  the  fact  that  the  wires  do  not  explode  the  same  on  every  shot.  The 
variation  in  the  geometry  of  the  initial  plasma  produced  by  the  exploding  wires  may  lead 
to  differences  in  the  quality  of  the  pinches. 
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4.2.  Pinhole  Camera  Results. 


The  pinhole  pictures  obtained  at  the  Phoenix  PRS  are  shown  in  Figure  4-2.  The  image  on 
the  left  is  a  digitized  scan  of  the  film.  Image  processing  was  used  to  enhance  the  contrast 
of  these  pictures  to  try  and  reveal  more  structure.  The  two  pictures  to  the  right  of  the 
original  scanned  image  have  increasing  contrast  going  to  the  right.  The  difference  in 
appearance  between  shot  numbers  1 1 1 8  -  1 125  and  numbers  1 126  -  1 129  is  due  to  the 
difference  in  the  pinholes  used  and  the  distance  of  the  film  from  the  pinhole  as  discussed 
in  Section  2.2.  To  get  more  accurate  information  about  structures  in  the  pinches,  an  EDP 
densitometer  was  used  to  obtain  contours  of  the  images. 

For  shot  numbers  1 1 18  -  1 125  the  densitometer  was  used  with  a  6.7  pm  aperture  and  an 
EDP  magnification  of  20.  Figure  4-3  shows  the  resulting  contour  image  for  shot  number 
1119.  From  this  scan,  the  length  of  the  hottest  part  of  the  plasma  from  the  pinhole  picture 
corresponds  to  1 .5  cm.  The  size  of  the  most  dense  structure  in  the  pinch  is  approximately 
350  -  500  pm,  but  since  this  corresponds  to  the  size  of  the  pinhole  used  of 400  pm,  it  is 
clear  that  the  picture  is  limited  by  the  resolution  of  the  camera.  In  anticipation  of  these 
results,  shot  numbers  1 126  -  1129  were  taken  with  a  smaller  pinholes  of  diameter  25  pm. 
For  the  scans  of  these  pictures  a  3.4  pm  aperture  was  used  with  an  EDP  magnification  of 
20.  The  resulting  contour  for  shot  number  1 126  is  shown  in  Figure  4-3.  The  length  of 
the  hottest  part  of  the  plasma  from  this  pinhole  picture  corresponds  to  1 .8  cm.  These 
scans  do  not  show  the  structure  that  was  apparent  in  the  scans  for  shot  numbers  1118- 
1125.  Since  film  saturation  was  suspected  to  be  a  possible  explanation  for  the  lack  of 
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Figure  4-2.  Pinhole  camera  pictures. 
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(a)  Shot  number  1119 


(b)  Shot  number  1126 


Figure  4-3.  EDP  scans  of  pinhole  pictures. 
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structure,  the  opacity  of  the  film  was  checked.  The  maximum  opacity  of  the  film  was 
found  to  be  2.2  which  corresponds  to  a  region  of  saturation  on  the  film  (Appendix  E).  It 
was  therefore  not  possible  to  use  this  data  to  determine  the  size  of  the  “hot  spots”.  The 
best  that  we  could  get  from  the  pinhole  data  was  that  the  size  of  the  “hot  spots”  were 
smaller  than  400  pm. 

4.3  Four  Frame  Camera  Results. 

The  pictures  obtained  using  the  four  frame  camera  are  shown  in  Figure  4-4.  Scans  of  the 
pictures  were  taken  using  the  EDP  densitometer.  Figure  4-5  shows  a  sample  of  the 
results  for  frame  3  of  shot  1 127  at  50x  EDP  magnification.  The  size  of  the  most  dense 
region  for  this  shot  is  approximately  500  pm,  corresponding  to  the  size  of  the  pinhole. 
These  results  verify  that  a  400  pm  pinhole  is  too  large  to  resolve  the  structure  of  the 
pinches.  Insufficient  resolution  was  found  to  be  a  problem  in  all  of  the  shots.  For  shot 
number  1 129  a  100  pm  pinhole  was  used  to  try  and  improve  the  results.  The  image 
appears  to  have  lower  resolution  than  the  pictures  obtained  with  the  400  pm  pinholes. 
Since  only  one  shot  with  this  size  pinhole  was  obtained  it  is  difficult  to  interpret  the  data. 
The  image  could  be  poor  because  of  the  camera  resolution,  or  the  pinch  could  have  been 
caught  in  expansion.  These  results  suggest  that  in  order  to  obtain  better  results  it  is 
necessary  to  use  smaller  pinholes.  The  use  of  smaller  pinholes  must  be  accompanied  by 
moving  the  camera  closer  to  the  plasma  so  that  the  total  radiation  to  the  film  can  be 
maintained.  The  resolution  of  the  four  frame  camera  pictures  do  not  give  sufficient 
resolution  to  determine  the  final  size  of  the  pinch. 
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4.4  Plasma  Parameters  from  Johann  Spectrometer  Profiles. 

The  spectra  used  in  this  analysis  were  scanned  using  a  Joyce  Loebl  scanning  densitometer 
(Appendix  F).  Between  4  and  6  scans  across  the  central  region  of  the  spectra  were  made 
and  averaged  with  0.5  mm  spacing  between  scans.  This  was  done  to  reduce  the  noise  due 
to  the  graininess  of  the  film.  The  Biomax  MR  film  used  for  the  Johann  and  flat  crystal 
spectrograph  was  calibrated  to  find  the  region  of  linear  response  (Appendix  E)  to  ensure 
that  film  saturation  was  not  a  problem. 

All  of  the  spectra  obtained  from  the  Phoenix  source  were  examined  to  determine  which 
shots  represented  good  useable  data.  When  the  spectrometers  are  aligned  properly  the 
spectra  should  appear  centered  on  the  film  both  horizontally  and  vertically.  Horizontal 
displacement  indicates  that  the  film  is  rotated  off  of  the  Rowland  circle  of  the  crystal, 
resulting  in  defocusing.  Vertical  displacement  on  the  film  indicates  that  the  x  rays  are  in 
a  plane  which  is  tilted  with  respect  to  the  plane  of  the  spectrometer.  Since  the  1  m 
diameter  of  the  Rowland  circle  used  is  large,  vertical  displacements  of  a  few  centimeters 
on  the  film  could  not  lead  to  significant  differences  in  the  geometry.  The  aberrations 
produced  by  this  displacement  are  discussed  by  Jones,  Paschen,  and  Nicholson.(30)  A 
secondary  problem  with  the  displacement  of  spectra  on  the  film  was  that  if  the  film  was 
too  far  displaced,  the  middle  of  the  line  did  not  appear  on  the  film.  Spectra  were  chosen 
which  were  centered  both  horizontally  and  vertically  within  1  cm  on  the  film.  Figures 
4-6  and  4-7  shows  the  location  of  the  extremes  of  the  spectra  on  the  film  horizontally  and 
vertically. 
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Figure  4-6.  Horizontal  position  of  spectra  on  film. 
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Eleven  “good”  spectra  were  obtained  for  the  Lp  line.  For  the  parallel  configuration  the 
shot  numbers  were  1 1 19  -  1122  and  for  the  perpendicular  configuration  the  shot  numbers 
were  1116,  and  1 1 18  - 1123.  Three  “good”  spectra  were  also  obtained  for  Ly.  One  was 
obtained  on  shot  1 124  for  the  parallel  configuration  and  two  were  obtained  on  shots  1124 
and  1125  for  the  perpendicular  configuration. 

Figure  4-8  (a)  and  (b)  show  the  typical  spectra  obtained  for  the  Lp  and  Ly  lines  on  the 
Phoenix  PRS.  The  Lp  spectra  has  two  A1 XII  lines  on  the  red  wing  of  the  line.  For  this 
reason  the  blue  wing  of  Lp  was  used  for  analysis.  The  spectra  for  Ly  shows  Ly,  Ls  and  LE 
in  addition  to  a  Si  XIII  line.  The  silicon  line  resulted  from  the  use  of  silicon  oil  to  clean 
the  machine.  Since  these  lines  occur  on  the  blue  side  of  the  line,  the  red  side  of  the  Ly 
line  was  used  for  analysis. 

Figure  4-9  shows  the  blue  side  of  the  Lp  lines  and  Figure  4-10  shows  the  red  wings  of  the 
Ly  lines  for  the  data  that  was  analyzed.  The  data  for  the  shots,  where  two  states  of 
polarization  were  obtained,  show  differences  in  half  widths  between  the  two  states  of 
polarization,  with  the  exception  of  shot  number  1 124.  Since  the  center  of  the  line  profiles 
are  strongly  determined  by  density  and  less  by  the  oscillating  electric  fields,  the 
differences  in  the  line  profiles  indicate  that  the  spectrometers  were  looking  at  regions 
with  different  densities.  The  fact  that  the  spectrometers  might  have  seen  different 
densities  in  the  plasma  is  reasonable  if  we  consider  the  difference  in  geometry  of  the 
spectrographs.  For  the  spectrograph  in  the  parallel  configuration,  a  scan  perpendicular 
tothe  dispersion  on  the  film  corresponds  to  viewing  different  regions  of  the  plasma  along 
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Figure  4-8.  Spectra  obtained  on  the  Phoenix  PRS. 
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Figure  4-10.  Red  wing  of  L,  lines. 
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the  length.  For  the  perpendicular  spectrograph,  a  scan  perpendicular  to  the  dispersion 
corresponds  to  a  scan  across  the  width  of  the  plasma.  Due  to  the  nature  of  the  elongated 
source,  at  best  one  can  expect  to  view  regions  with  similar  average  densities. 

Once  the  spectra  for  the  blue  wing  of  Lp  and  the  red  wing  of  Ly  was  obtained  for  the 
“good”  shots,  the  model  parameters  that  produced  these  lines  could  be  found.  Since  there 
are  six  parameters  in  the  model  this  corresponds  to  a  six  dimensional  search.  Each 
parameter  can  vary  over  a  range  of  approximately  5,  so  it  would  take  approximately  1010 
searches  to  cover  the  entire  range  to  a  resolution  of  0. 1 .  To  avoid  having  to  search  the 
entire  range,  a  six  dimensional  Sobol  quasi-random  sequence  was  employed. (32)  A  Sobol 
sequence  generates  numbers  in  an  n  dimensional  space  such  that  each  additional  point 
added  is  further  in  the  space  from  previous  points.  As  more  of  the  sequence  is  produced, 
the  space  is  searched  to  a  finer  scale.  This  is  best  illustrated  by  looking  at  the  numbers 
generated  for  a  two  dimensional  Sobol  sequence.  Figure  4-11  shows  the  two  dimensional 
Sobol  sequence  up  to  1024  points. 

A  program  was  written  to  find  the  parameters  which  matched  the  experimental  lines. 
Eleven  points  were  taken  from  each  of  the  experimental  profiles  to  be  compared  to  the 
profiles  generated  by  the  model.  The  scheme  was  to  take  a  set  of  parameters  generated 
by  the  Sobol  sequence,  produce  the  profile  at  these  points,  and  then  compare  them  to  all 
sets  of  experimental  data  for  the  line.  If  all  1 1  points  in  the  experimental  profile  matched, 
to  within  a  set  tolerance,  the  parameters  were  recorded  as  a  solution  for  the  experimental 
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data.  This  was  done  for  a  number  of  possible  baselines  to  ensure  that  the  results  were  not 
a  function  of  a  particular  baseline  choice. 

For  the  Lp  and  Ly  lines,  1 1  different  baselines  were  tested.  In  Figure  4-9  and  4-10  the 
baselines  were  estimated  from  where  the  wings  of  the  lines  died  out.  This  represents  a 
best  initial  guess  of  the  baselines.  The  height  of  these  lines  were  normalized  to  unity  and 
baselines  were  adjusted  in  terms  of  a  percentage  of  this  height.  Raising  the  baseline  by 
10%,  for  example,  corresponded  to  subtracting  0.1  from  the  best  guess  profile  and 
renormalizing  the  line.  Baselines  ranging  from  raising  the  baseline  20%  to  lowering  the 
baseline  by  30%  were  tested  by  increments  of  2.5%  for  the  Lp  and  Ly  lines. 

The  range  of  model  parameters  tested  ranged  from  2.0  x  104  to  7.0.  For  Lp,  25,000 
searches  were  performed  and  for  Ly  100,000  searches  were  performed.  The  number  of 
searches  was  based  on  the  number  of  solutions  that  were  produced.  Figure  4-12  shows  a 
sample  of  the  results  from  the  search  of  the  Lp  line  for  shot  number  1 1 1 8  in  the 
perpendicular  configuration  and  Figure  4-13  shows  the  results  from  the  search  of  the  L., 
line  for  shot  number  1 124  in  the  parallel  configuration.  Figure  4- 14(a)  and  4- 14(b)  have 
the  fits  of  the  produced  profiles  to  the  experimental  profiles  for  these  shots.  These 
searches  illustrate  that  there  are  many  sets  of  parameters  which  can  reproduce  the 
experimental  data.  A  list  of  the  range  of  each  parameter,  where  solutions  have  been 
found,  is  summarized  in  Table  4-2  for  the  Lp  line  and  Table  4-3  for  the  Ly  line. 


82 


83 


I  1  i  T  I  1  i  1  i  1  i  1  i  1  i 

"'O'n^renc^  —  3 

snpjA  isisurejEj 


84 


Table  4-2.  Range  of  parameters  for  Lyman  (3  which  fit  the  experimental  profiles. 


b 

c 

8 

£ 

z 

t 

16P 

1.1  -2.8 

.02 -.18 

1.4 -3.7 

1.4 -5.3 

.16-5.4 

1.3 -2.2 

18P 

0.9  -  4.4 

.02  -.26 

0.5  -.44 

1.9 -6.6 

0.5  -  5.3 

0.7 -4.4 

19L 

2.2 -3.9 

.02  -.15 

1.1  -4.2 

1.2 -6.9 

0.2  -  5.0 

.02-1.6 

19P 

1.2 -3.3 

.02  -.11 

1.2 -4.4 

2.4  -  6.6 

0.2 -5.3 

1.5 -3.8 

20L 

0.9 -3.1 

.02  -.33 

1.1  -3.3 

1.2 -6.0 

0.8  -  5.0 

0.2  -  3.4 

20P 

1.0 -2.9 

.03  -.22 

1.2 -4.4 

1.8 -6.8 

0.2  -  5.4 

1.8 -3.9 

21L 

3.0 -3.5 

.02  -.04 

2.2 -4.1 

4.2  -  6.4 

0.7 -4.7 

0.2  -  0.8 

21P 

1.3  -3.3 

.02  -.15 

0.8 -4.2 

2.7  -  6.8 

0.6  -  5.3 

1.1  -2.8 

22L 

1.0 -3.6 

.004  -.49 

0.5 -4.5 

1.1  -6.2 

0.6 -5.1 

.03  -  3.0 

22P 

2.1  -3.3 

.03  -.33 

0.9  -  3.3 

0.9  -  5.6 

0.1  -4.9 

0.6 -3.0 

23P 

1.7 -3.9 

.01  -.47 

0.9  -  4.5 

1.3 -6.9 

0.1  -5.4 

0.3 -3.1 

Table  4-3.  Range  of  parameters  for  Lyman  y  which  fit  the  experimental  profiles. 


b 

c 

g 

£ 

z 

t 

24L 

1.5  -2.9 

.01  -.5 

1.0 -4.3 

0.5 -6.1 

0.1  -  1.9 

0.1  -4.0 

24P 

1.3 -2.2 

.01  -0.2 

0.4  -  4.2 

0.6  -  6.7 

.05  -  2.2 

0.2  -  4.4 

25P 

1.3  -2.4 

.05  -.3 

1.0 -4.3 

0.9 -3.2 

.03  -.22 

0.6  -  4.07 

It  is  clear  from  these  results  that  a  fit  to  one  line  does  not  produce  a  unique  solution. 
Since  there  are  effectively  seven  degrees  of  freedom,  six  parameters  and  one  choice  of 
baseline,  less  parameters  or  more  constraints  are  needed  to  find  unique  parameters.  In 
order  to  get  an  estimate  of  the  key  parameters  from  the  Lp  and  Ly  lines,  the  searches  were 
performed  again  using  fewer  of  the  parameters. 
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A  reduced  search  for  the  Lp  lines,  with  electric  field  anisotropy  turned  off,  was  performed 
to  get  an  estimate  of  the  plasma  parameters  b,  C,  g,  e,  and  t.  We  expect  the  lines  to 
depend  strongly  on  s,  since  the  lines  appear  to  be  predominantly  broadened  by  the 
oscillating  electric  fields.  This  is  illustrated  by  looking  at  the  plots  of  log  intensity  vs.  AX 
for  the  Lp  lines  (Figure  4-15).  These  plots  show  a  quasi-linear  behavior  indicative  of  a 
superposition  of  Gaussian  profiles.  The  superposition  of  Gaussian  profiles  results  from 
the  sum  of  the  envelopes  of  the  Stark  broadened  components  under  the  action  of 
oscillating  electric  fields.  The  search  performed  with  the  reduced  parameter  set  for  all  Lp 
lines  lead  to  the  parameters  b  =  2.0  -  4.0,  C  =  0.02  -  0.4,  g  =  0.3  -  2.5,  s  =  1.0  -  2.5,  and 
t=  1-3. 

A  reduced  search  for  the  gamma  line  was  also  performed.  Since  the  Ly  line  has  a  central 
component  which  is  unaffected  by  Stark  broadening,  the  center  of  the  line  is  dominated 
by  Doppler  broadening.  The  oscillating  electric  field  will  dominate  in  the  wing  of  the 
line  where  the  Doppler  profile  of  the  central  component  dies  out.  For  this  reason,  the 
wings  of  the  gamma  line  were  fitted  with  the  reduced  parameter  set  C,  g,  s,  and  z.  The 
parameter  z  was  included  to  determine  the  anisotropy  of  the  electric  fields  from  shot 
number  1124  where  the  spectrometers  are  looking  at  a  similar  region  of  the  plasma. 

Several  sets  of  the  reduced  parameters  were  found  to  fit  the  wings  of  the  lines.  To  reduce 
the  number  of  solutions,  we  required  that  there  be  a  reciprocal  relationship  for  z  in  the 
perpendicular  configuration  to  the  parallel  configuration  for  the  spectrometers,  since  the 
meaning  of  a  and  b  in  z  are  swapped  in  the  two  configurations.  The  results  for  the 
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reduced  Ly  search  of  the  wings  were:  C  =  0.02  -  0.20,  g  =  0.3  -  2.5,  s  =  1  -  2.5,  z  =  1  -  5. 

These  results  show  agreement  with  the  results  obtained  from  the  reduced  search  of  the  Lp 
lines. 

The  additional  constraints  needed  in  the  full  model  can  be  provided  by  having  the  model 
fit  multiple  lines  with  the  same  parameters.  The  flat  crystal  spectrograph  provided  the 
data  for  which  this  was  done. 

4.5  Plasma  Parameters  from  Flat  Crystal  Spectrometer  Profiles. 

Figure  4-16  shows  the  spectra  obtained  for  the  flat  crystal  spectrograph  at  the  Phoenix 
PRS.  The  search  for  the  model  parameters  which  would  reproduce  the  lines  was  carried 
out  in  the  same  manner  as  for  the  Johann  spectrometer  data.  The  additional  constraint  in 
the  search  was  that  the  same  parameters  had  to  produce  all  of  the  lines  in  the  spectra. 
However,  opacity  was  only  applied  to  the  L?  line,  since  opacity  decreases  rapidly  along 
the  spectral  series.  The  central  parts  of  the  lines  were  used  for  the  search,  since  there  is 
an  overlap  in  the  wings  of  the  lines  of  the  flat  crystal  spectra.  By  fitting  only  the  central 
part  of  the  lines,  many  solutions  were  found.  The  criteria  used  to  determine  if  the  found 
profiles  were  acceptable  was  that  the  region  of  the  overlapping  lines  between  L§  and  L£ 
had  to  match  the  experimental  profile.  The  fits  which  contained  the  three  lines  Ly,  L5  and 
Le  were  the  only  ones  which  led  to  a  unique  set  of  parameters  and  this  resulted  only  for 
the  baseline  raised  5%  from  the  best  guess  baselines  given  in  Figure  4-16.  The  two 
spectra  which  contained  these  three  lines  were  shot  number  1116  and  shot  number  1 120. 
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A  best  fit  for  the  model  parameters  in  shot  number  1116  was  found  to  be  b  =  1.5  -  2.5,  C 
=  0.01  -0.4,  g  =  0.5- 1.5,  £  =  0.1  -2.5,2=  1-4.5,  and  t=  0.7  - 1.8  and  a  best  fit  for  the 
model  parameters  in  shot  number  1 120  was  found  to  be  b  =  0.9  -  2,  C  =  0.01  -  0.2,  g  =  0.7 
-  2,  £  =  0.05  - 1.7,  z  =  1-  5,  and  t  =  1.5  -  4.5.  The  experimental  spectra  for  shot  number 
1116  and  1 120  are  shown  in  Figure  4-17  (a)  and  (b)  with  one  set  of  parameters  which  fit 
the  spectra.  From  the  fits  of  the  Lp,  Ly  and  flat  crystal  profiles  the  parameters  in  the 
Phoenix  source  were  consolidated. 

4.6  Interpretation  of  the  Plasma  Parameters  Found. 

The  parameter  b  was  found  to  be  approximately  1 .0  -  2.5.  This  corresponds  to  a  bulk 
plasma  velocity  of  3.0  -  7.5  x  107  cm  /  s.  This  result  is  reasonable  since  velocities  of  this 
order  have  been  observed  on  many  other  similar  z-pinch  devices.(8) 

The  parameter  C  was  found  to  have  values  from  0.01  to  0.40.  This  corresponds  to  a 
reasonable  electron  density  of  0.3  -  6  x  1021  cm'3. 

The  parameter  g  =  0.5  to  2.5  corresponds  to  a  power  spectrum  with  a  width  ranging  from 
0.25  <X)pe  to  1.25  c ope.  Values  of  g  higher  than  2.5  were  excluded,  since  this  corresponds  to 
a  width  greater  than  the  plasma  frequency,  and  is  not  physical. 

The  value  s  obtained  from  the  Johann  spectrometers  for  the  Lp  and  Ly  profiles  was  1  -  2.5 
and  the  value  of  s  obtained  from  the  flat  crystal  spectrograph  for  the  Ly,  L6,  and  Le 
profiles  was  0.05  -  2.5.  The  actual  value  of  s  is  probably  somewhere  between  these  two 
results  in  the  range  1  -  2,  since  there  is  some  uncertainty  associated  with  each  of  these 
two  results.  The  uncertainty  in  the  results  from  the  Johann  spectrometers  comes  from 
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5.40  5.45  5.50  5.55  5.60  5.65  5.70  5.75  5.80  5.85  5.90 


using  a  reduced  parameter  set  to  obtain  an  estimate  of  s,  while  the  uncertainty  in  the 
results  from  the  flat  crystal  spectrometer  comes  from  the  inherent  low  resolution  of  the 
spectrometer.  The  values  of  £  =  1  -  2  and  C  =  .01  -  0.4  correspond  to  an  oscillatory 
electric  field  of  E0»  1-9  GV  /  cm.  These  values  for  the  oscillating  electric  fields  are 

reasonable  in  terms  of  energy  considerations.  The  energy  density  in  the  electric  fields  of 
E02  /  16tc  was  compared  to  the  thermal  energy  density  NeT  and  the  kinetic  energy  density 
KE  =  (3 N'.vJ)  /  2  in  the  plasma.  The  energy  density  of  the  oscillations  were  found  to  be 
approximately  0.5-4  times  that  of  the  thermal  energy  density,  by  assuming  an  electron 
temperature  from  500  keV  to  1000  keV.  The  high  end  values  of  this  ratio  seem 
physically  improbable.  Because  the  waves  have  an  energy  density  comparable  to  the 
thermal  energy  density,  there  is  a  strong  Langmuir  turbulence.  The  ratio  of  the  energy 
density  of  the  oscillations  to  the  kinetic  energy  density  varies  from  0.01  to  0.3,  indicating 
that  the  energy  density  balance  was  not  violated. 

The  parameter  z  =  1  -  5  corresponds  to  an  anisotropy  of  the  electric  field  with  the  greatest 
field  strength  along  the  direction  of  the  bulk  current.  This  agrees  with  the  theory  which 
suggests  that  these  anomalous  electric  fields  are  produced  by  runaway  electrons  in  the 
discharge.  The  runaway  electrons  are  expected  to  be  predominantly  along  the  direction 
of  the  bulk  current. (34) 

The  parameter  t  showed  a  variation  from  1  and  3  for  the  Lp  line  and  0.5  to  2  in  the  flat 
crystal  results  for  the  Ly  line.  These  are  reasonable  optical  depths  for  the  Lp  and  Ly  lines. 
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SECTION  5 

DISCUSIONS  AND  CONCLUSIONS 


By  using  several  spectrometers  we  were  able  to  measure  the  level  and  angular 
distribution  of  oscillating  electric  fields  in  a  superdense  plasma  produced  on  the  Phoenix 
device.  Since  the  energy  density  in  these  fields  is  on  the  order  of  the  thermal  energy 
density,  this  suggests  that  these  oscillations  may  lead  to  the  instability  which  prevents  the 
highest  ionization  states  of  higher  Z  elements  from  being  reached  in  plasma  radiation 
sources.  Further  experiments,  using  more  sophisticated  designs  are  necessary  to 
determine  the  field  values  more  accurately. 

There  are  several  improvements  which  can  be  made  to  measure  the  two  states  of 
polarization  of  a  line  profile.  One  of  the  problems  encountered  on  the  Phoenix 
experiment  was  that  the  two  Johann  spectrometers  were  not  guaranteed  to  be  looking  at 
the  same  part  of  the  plasma.  This  resulted  from  the  extension  of  the  source  along  the 
wires  and  the  placement  of  the  spectrometers  on  opposite  sides  of  the  source. 

The  difficulty  of  having  an  extended  source  is  that  the  resulting  spectra  is  an  average  over 
the  length  of  the  source.  A  point  source  is  necessary  to  measure  spectra  in  two  states  of 
polarization  from  a  localized  region  of  the  plasma. 

The  problem  associated  with  the  location  of  the  spectrometers  could  be  resolved  if  both 
states  of  polarization  could  be  extracted  from  the  same  line  of  sight  to  the  plasma. 
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For  a  point  source,  such  a  spectrometer  can  be  designed  by  using  a  thin  flat  crystal  as  a 
kind  of  beam  splitter  as  shown  in  Figure  5-1.  In  the  ideal  case,  if  a  crystal  which  leads  to 
a  Bragg  angle  of  45°  can  be  found,  for  the  desired  wavelength,  and  can  be  made  thin 
enough  so  that  enough  radiation  can  pass  through,  then  a  spectrometer  which  completely 
polarizes  the  x  rays  can  be  made.  Figure  C-3.  shows  the  calculation  of  transmission  and 
reflection  of  6.053  A,  for  two  states  of  polarization  at  45.3°,  for  thin  crystals  of  N  layers 
of  quartz  10K).  These  curves  show  that  at  45°  the  n  polarization  is  not  reflected  with  any 
significant  intensity.  More  importantly,  the  n  component  is  transmitted  without  any 
difference  than  that  for  bulk  material  absorption.  This  means  if  we  set  the  thin  crystal  A 
close  to  45°,  such  that  only  the  a  component  is  reflected,  then  we  can  use  a  second  crystal 
B,  also  oriented  at  45°,  but  in  the  opposite  state  of  polarization,  to  Bragg  reflect  the  n 
component  that  passed  through  the  thin  crystal  A.  In  doing  so,  we  then  have  extracted 
both  components  of  polarization  from  the  same  line  of  sight.  In  practice  however,  it  will 
not  always  be  possible  to  find  a  crystal  with  the  correct  2d  value  which  can  be  cleaved 
thin  enough  to  transmit  sufficient  radiation  at  the  desired  wavelength,  especially  for 
longer  wavelengths  where  absorption  coefficients  are  large. 

Even  if  we  cannot  find  a  thin  crystal  A,  which  results  in  a  Bragg  angle  at  45°,  we  can  still 
measure  the  profiles  in  two  states  of  polarization  if  we  can  find  a  thin  enough  crystal  (A) 
where  the  Bragg  angle  results  in  a  difference  in  the  amount  of  radiation  reflected  for  the  n 
and  a  components.  From  calculations  of  thin  crystals,  we  know  the  amount  of  the  n  and 
cr  components  that  will  be  reflected  from  the  thin  crystal  into  the  first  detector.  The 
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radiation  passing  through  the  crystal  will  still  be  unpolarized  for  the  most  part.  The 
second  crystal  B  must  still  have  a  Bragg  angle  close  to  45°  so  that  only  one  state  of 
polarization  makes  it  to  the  second  detector.  By  knowing  the  relative  ratio  of  reflected 
intensities  for  crystal  A  and  having  a  spectra  in  one  polarization,  the  a  and  n  spectra  can 
be  extracted. 

In  addition  to  improving  the  geometry  of  the  spectrometers,  it  would  also  be  desirable  to 
have  time  resolved  spectra.  This  could  be  achieved  by  using  x-ray  streak  cameras  for  the 
detectors.  The  one  disadvantage,  however,  of  using  streak  cameras  is  the  reduction  of 
detector  resolution.  Still,  spectral  profiles  at  various  stages  of  the  pinch  event  would 
provide  detailed  information  about  the  time  evolution  of  the  plasma. 

It  would  also  be  desirable  to  have  a  high  resolution  spectrometer  to  look  at  the  Lyman 
lines,  instead  of  the  flat  crystal  spectrometer  used  in  this  experiment.  A  Johann  type 
spectrometer  could  be  used  in  place  of  the  flat  crystal  spectrometer.  With  a  higher 
resolution  spectra  of  the  Lyman  series,  the  parameters  in  the  plasma  could  be  determined 
with  more  accuracy. 

In  addition  to  the  improvements  made  to  the  spectrometers,  several  improvements  can  be 
made  to  the  other  diagnostics  used  on  the  Phoenix  source.  Several  improvements  for 
measuring  the  size  of  pinch  can  be  made  on  the  pinhole  camera.  First,  it  is  necessary  to 
use  smaller  pinholes  in  the  camera.  This  would  bring  the  resolution  up  high  enough  so 
that  the  pinches  could  be  resolved.  In  addition  to  using  smaller  pinholes  it  is  also 
necessary  to  determine  the  optimum  filters  for  the  cameras.  As  we  saw  from  the  pinhole 
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pictures  on  the  Phoenix  experiment,  if  too  little  filtration  was  used,  the  pictures  became 
saturated.  Ideally  the  intensity  of  the  source  should  be  calibrated  using  known  filter 
thicknesses.  This  would  require  being  able  to  get  an  absolute  optical  density  from  a  few 
pictures  with  known  filters  in  place.  A  starting  point  for  changing  filters  and  pinhole 
sizes  would  then  be  established.  It  would  also  be  ideal  to  use  pinholes  made  out  of  a  high 
Z  material,  like  platinum,  to  prevent  radiation  from  passing  through  the  pinhole  substrate. 
These  improvements  would  lead  to  the  best  results  for  pinch  size  determination  using  a 
pinhole  camera. 

The  suggested  improvements  on  the  regular  pinhole  camera  can  also  be  made  to  the  four 
frame  camera.  In  addition  to  using  smaller  pinholes  and  the  appropriate  filtration,  the 
intensity  at  the  camera  has  to  be  maintained,  since  each  frame  of  the  camera  is  only  open 
for  several  nanoseconds.  The  intensity  at  the  detector  can  be  increased,  to  offset  the 
reduction  in  intensity  from  using  smaller  pinholes  by  moving  the  camera  closer  to  the 
source.  With  the  correct  pinholes,  filters,  and  intensity  at  the  detector,  many  shots  would 
then  have  to  be  made  to  get  enough  information  to  determine  the  size  of  the  pinches. 
Pictures  would  also  have  to  be  taken  at  several  places  in  the  radiation  pulse  to  ensure  the 
smallest  pinches  are  observed. 

With  these  improvements  to  the  diagnostics,  a  better  understanding  of  superdense 
plasmas  like  the  Phoenix  device  could  be  achieved  and  this  could  lead  to  improved 
performance  of  pulsed  radiation  sources. 
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APPENDIX  A 

CONVOLUTION  OF  LINE  PROFILES 


A.1  CONVOLUTION  OF  DOPPLER,  HOLTSMARK  AND  LORENTZ 
PROFILES. 

To  arrive  at  the  convolution  of  the  Doppler,  Holtsmark,  and  Lorentz  profiles  we  first 
convolute  the  Doppler  and  Holtsmark  profiles  and  then  convolute  the  result  with  the 
Lorentz  profile. 

A.1.1  Convolution  of  the  Gaussian  and  Holtsmark  Profiles. 

We  are  first  interested  in  convoluting  the  Doppler  profile  given  by 
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A.1.2  Convolution  of  the  Gauss-Holtsmark  and  Lorentz  Profiles 

Once  we  have  the  convolution  of  the  Gaussian  and  Holtsmark  profiles  we  need  to 
convolute  this  profile  with  the  Lorentz  profile  given  by 
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Making  the  substitution 
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the  integral  becomes 
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The  integrals  in  z  both  contain  two  poles.  We  must  choose  the  contours  correctly  such 
that  the  part  of  contour  off  of  the  real  axes  vanishes  as  z  goes  to  ±  oo.  For  y  >  0  we  have 
to  make  the  integration  with  the  imaginary  part  of  z  negative.  For  y  <  0  we  have  to  use  a 
positive  value  for  z.  Applying  the  residue  theorem,  the  integral  becomes 
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The  sign  of  the  first  term  is  positive  because  a  minus  sign  is  picked  up  by  the  direction  of 
the  contour.  Reversing  the  limits  of  integration  and  letting  y  become  -y  in  the  second 
integral  we  have 
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Putting  the  exponential  back  in  terms  of  sin  and  cosine  and  taking  the  real  part  we  arrive 
at  the  final  integral  for  the  profile 
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APPENDIX  B 

DISPERSION  OF  FLAT  CRYSTAL  SPECTROMETER 


B.1  CALCULATION  OF  DISPERSION  FOR  THE  AUBURN  FLAT 
CRYSTAL  SPECTROMETER 


To  calculate  the  dispersion  X(x)  on  the  film  for  the  flat  crystal  spectrometer,  we  need  to 
define  some  physical  dimensions  for  the  spectrometer.  The  setup  for  the  Auburn  flat 
crystal  spectrometer  is  shown  in  Figure  3-13.  We  take  the  distance  from  the  crystal  to  the 
source  to  be  U  and  the  distance  from  the  center  of  the  two  crystals  to  the  film  plane  to  be 
V.  The  flat  crystals  are  set  to  an  angle  a  relative  to  the  source.  With  these  parameters 
defined  we  can  use  geometry  to  find  the  dispersion. 

To  begin,  we  find  the  wavelength  X  of  the  source  that  will  satisfy  Bragg’s  law  a  distance 
h  from  the  center  of  the  two  crystals.  From  the  law  of  sines,  the  angle  (3  that  a  ray  makes 
a  distance  h  from  the  center  of  the  crystals  is  given  by 


P 


=  sin 


-i 


_ U  sin(q) _ 

yju2  +h2  +2hUcos(a ) 


(B.l) 


From  Bragg’s  Law,  X  =  2d  sin  [3 ,  we  find  the  wavelength  X  that  corresponds  to  this  angle 
is 
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Solving  for  h  we  get 
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We  take  the  positive  solution,  since  the  negative  root  corresponds  to  the  same  Bragg 
angle  but  on  the  opposite  side  of  the  normal  that  we  are  interested  in. 

We  can  now  proceed  to  find  x(k)  by  finding  the  intersection  of  the  Bragg  reflected  ray 
from  the  crystal  and  the  line  of  the  film.  The  equation  for  the  line  of  the  ray  is  given  by 


where  the  slope  m  is 


Yr  =  mXR  +  b, 
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m  =  tan(a  +  P) .  (B.5) 

The  y-intercept  b  can  be  found  since  we  know  that  the  point 

YR=  h  sin(a)  XR  =  h  cos(a)  (B.6) 

lies  on  the  reflected  line.  Solving  for  b  and  plugging  it  back  into  Eq.  (B.4)  we  have 

Yr  =  tan(a  +  P)^  +  /z[sin(a)  -  tan(a  +  p)  cos(a)] .  (B.7) 

This  line  intersects  the  film  when  YR  =  V .  Plugging  V  in  for  YR  in  Eq.  (B.3)  and  solving 
fox  XR  we  get 

XR=h  cos(a)  +  F~>’sin(a) .  (B.8) 

*  tan(a  +  p) 

x(A.)  can  then  be  found  by  putting  in  Eq.  (B.l)  for  h,  that  is 

x(X)  =  U  cos(a)%  +  F~^Sm(a)  (B.9) 
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To  find  A.(x) ,  the  expression  for  x(X)  is  most  easily  inverted  numerically. 
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APPENDIX  C 

POLARIZATION  BY  CRYSTALS 


C.l  POLARIZATION  CALCULATIONS  FOR  REFLECTION  FROM  SEMI¬ 
INFINITE  CRYSTALS. 

To  properly  design  spectrographs  that  are  capable  of  distinguishing  between  states  of 
polarization  it  is  necessary  to  understand  how  x  rays  are  polarized  under  Bragg  reflection 
from  crystals.  The  problem  of  x  ray  scattering  from  crystals  is  treated  in  detail  by  Henke 
et.  al  using  modified  Darwin  Prin  theory.(35)  The  calculations  based  on  the  MDP  theory 
will  be  summarized  and  applied  here  to  the  quartz  crystals  used  in  the  experiment. 

From  MDP  theory  the  reflected  intensity  for  a  semi-infinite  crystal  is  given  by 

|2 

J__  _ 

4  (a  +  e)  ± 

The  parameters  s  and  a  are  defined  such  that  the  fractional  reflected  amplitude  is  -is  and 
the  fractional  transmitted  amplitude  is  1-/ct  as  shown  in  Figure  C-l .  s  and  a  can  be 
calculated  by 
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where 
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is  the  number  of  unit  cells  per  unit  area  M  multiplied  by  the  average  structure  factor  F  of 
the  unit  cell.  The  average  structure  factor  F  is  calculated  by  averaging  the  contributions 
of  the  x  rays  scattered  from  each  atom  in  the  unit  cell.  The  quantity  mq  is  the  number  of 
atoms  of  type  q  that  are  a  distance  zq  from  the  reference  plane  of  scattering  as  per  Figure 
C- 1 .  The  scattering  from  each  type  of  atom  is  given  by  an  atomic  scattering  factor 

/  =  -A(0)+*/2(0)-A/(G),  (C.5) 


which  has  an  angle  independent  part  described  by  scattering  in  the  forward  direction  and 
an  angle  dependent  correction  term  Af(0)  due  to  the  interference  affects  of  waves 
scattered  from  different  parts  of  the  atom.  The  angle  independent  scattering  is  primarily 
due  to  photoabsorption  and  Compton  scattering.  The  angle  dependent  correction  term  is 
given  by 


Af(Q)=Z-f0  (C6) 

wher ef0  is  the  atomic  scattering  factor  which  depends  on  sinQ  /  X.  The  atomic  scattering 
factors  for  atoms  can  be  found  in  tables.(36) 

The  factor  s  ,  given  by 


27idsin0 

s  = - rrm , 

^  9 


( C.7 ) 


represents  the  phase  difference  from  the  phase  required  for  the  nF  order  reflection.  Order 
m  =  0  corresponds  to  non  Bragg  reflection,  m  =  1  to  1st  order  Bragg  reflection,  etc. 
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The  polarization  factor  P(2Q)  for  the  incident  radiation  is  given  by  1  for  the  c  polarization 
and  cos(2Q)  for  the  n  polarization,  a  is  the  polarization  with  the  electric  field 
perpendicular  to  the  plane  of  incidence  and  n  is  the  component  with  polarization  parallel 
to  the  plane  of  incidence. 

A  factor  k  is  introduced  to  take  into  account  the  quality  of  the  crystal.  In  many  crystals 
there  can  be  a  misorientation  of  “crystal  segments”  as  illustrated  in  Figure  C-l .  A  value 
of  k  =  0  corresponds  to  a  mosaic  crystal  and  a  value  of  k  =  1  corresponds  to  a  perfect 
crystal.  Since  quartz  is  nearly  a  perfect  crystal,  in  all  successive  calculations  involving 
quartz,  k  will  be  taken  to  be  unity.(37) 

The  use  of  Modified  Darwin  Prin  theory  for  calculating  reflectivities  is  applicable  only 
for  0  >  0C  where 


ec*3X,1 
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When  the  angle  is  smaller  than  this  critical  angle,  the  Optical  Electromagnetic  Model 
must  be  used.  This  consists  of  defining  a  complex  index  of  refraction  in  terms  of  s  and  a 
and  solving  the  Fresnel  equations. 

C.1.1  POLARIZATION  OF  X  RAYS  FROM  QUARTZ  1010  IN  BRAGG 
REFLECTION. 


To  show  a  sample  of  how  MDP  theory  was  used  to  find  the  relative  reflectivity  of  the  a 
and  n  components  for  the  crystals  used  in  the  experiment,  we  calculate  the  reflection  of 
6.053 A  from  the  10J_0  plane  of  quartz  (Si02)3  under  Bragg  reflection. 


C-4 


Quartz  has  a  hexagonal  structure  with  the  dimensions  of  the  unit  cell  given  by 

a0  =  4.913  A  c  =  5.4  A.  (C.9) 

The  2d  spacing  in  quartz  for  the  10J.0  plane  is  8.51  A  and  the  critical  cutoff  angle  for 
using  MDP  theory  is  0C  =  0.86° .  The  position  of  the  atoms  in  the  unit  cell  of  quartz 
given  by  Wychoff  are(38) 

Si:  -0.465,  -0.465,  0.333  0.465,  0.000,  0.000  0.000,  0.465,  0.670 

O:  0.415,  0.272,  0.120  -0.143,  0.415,  0.455  -0.272,  0.143,  0.787 

0.143,-0.272,-0.120  0.272,0.415,0.0055  0.415,-0.143,0.215 

The  reflectivity  for  a  crystal  is  obtained  from  integrating  the  rocking  curve  given  by  Eq. 
(C.  1 ).  The  rocking  curve  is  the  reflected  intensity  I,  as  a  function  of  the  angle  6,  for  a 
given  wavelength.  In  order  to  calculate  the  rocking  curve  for  quartz,  we  have  to  first 
compute  MF(0)  and  MF(Q)  for  quartz  using  Eq.  (C.4). 

In  general,  the  distance  zq  of  an  atom  in  the  unit  cell  from  the  plane  of  scattering  is 
calculated  by  its  miller  indices  h,  k,  and  /  using  zq  /  d  =  hxq  +kyq  +  lzq ,  but  since  the 
1 0 1_0  plane  of  quartz  corresponds  to  planes  perpendicular  to  the  x  axes,  zq  in  the 
expression  for  MF  is  simply  given  by  the  x  position  of  the  atom  in  the  unit  cell.  A  further 
simplification  arises  in  calculating  MF  for  quartz  because  the  atoms  in  the  unit  cell  are 
symmetric  about  x  =  0.  This  allows  us  to  simplify  the  expression  for  MF  to 


(C.10) 


+  2kl(a2 bc){cos$ cosy  -  cosa)  +  2//z(aZ>2c)(cosy  cosa  -cosP)  ,(40)  (C.14) 

a,  b  and  c  are  points  where  the  scattering  plane  intersects  the  axis  and  a,  p  and  y  are  the 
direction  cosines  between  the  axes  in  a  particular  crystal  “system”.  For  quartz  which  is  a 
hexagonal  “system”,  a  =  P  =  90 ,  and  y  =  120 .  To  calculate  M  we  need  to  know  the 
Cartesian  dimensions  of  the  unit  cell  a,  b,  and  c.  For  hexagonal  crystals  a  =  b  =  a0  and 
we  find  that  M  =  3.766  x  1 014  cm'2. 
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To  finish  calculating  MF(0)  and  MF(Q)  we  need  the  scattering  factor  for  Silicon^,- and 
the  scattering  factor  for  Oxygen/,.  The  angle  independent  scattering  terms//  and/  in 
Eq.  (C.5)  are  found  in  tables  by  wavelength  for  each  element.(35)  For  the  angle  dependent 
term  A/(0)  we  need  to  calculate  sinQ  /  X  to  find  the  atomic  scattering  factors/  for  Silicon 


and  Oxygen.  For  quartz  1 010  at  the  wavelength  l  =  6.053  A  the  ratio  is  given  by 


sinG 


2d 


0.12. 


(C.15) 


The  atomic  scattering  factors  for  Silicon  and  Oxygen  and  were  found  in  the  tables  to  be 

A f(®)si =  U-35  and  Af(B)0  =  7.25 .(36)  The  scattering  factors  at  X  =  6.053  A  are  then  given 
by 


fsi=  12.33  +  3.58  /  -  (14-1 1.35)  =  9.68  +  3.58  i  (C.16) 

and 


f0=  8.31  +  0.51  i  -  (8-7.25)  =  7.56  +  0.51  i,  (C.17) 

resulting  in 

MF(0)  =  3.2x10 16  +5.2x10 15  /  (C.18) 

and 


MF(  0)  =  (-6.7x1 0 1 5  )  +  (-1.4x1 0 1 5  )i .  (C.  1 9) 

A  computer  program  was  used  to  compute  the  rocking  curves  from  Eq.  (C.l)  for  angles 
around  the  Bragg  angle  for  both  the  a  and  k  states  of  polarization.  These  curves  were 
then  integrated  to  yield  the  reflected  intensities  Ra  and  i?,.  The  ratio  of  the  reflected 
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intensities  for  the  n  and  a  polarization  gives  a  measure  of  how  the  crystal  polarizes 
incident  x  rays. 

The  resulting  rocking  curves  for  X  =  6.053  A  along  with  the  results  for  two  other 

wavelengths,  1 .50  A  and  4.25  A,  are  given  in  Figure  C-2.  The  ratio  of  reflected 

intensities  for  the  two  states  of  polarization  at  these  wavelengths  are  listed  in  table  IV. 

Table  C-l .  Ratio  of  integrated  reflectivity  for  c  and  n  polarizations  under  Bragg  reflection 
from  quartz  10J_0  for  several  wavelengths. 


Wavelength  (A) 

Angle  (degrees) 

K'K 

1.5 

10 

0.97 

4.255 

30 

0.33 

6.053 

45.2 

0.001 

Most  important  for  this  experiment  is  the  result  for  0  =  45°.  Here  only  one  state  of 
polarization  is  reflected  by  the  crystal. 

C.2  TRANSMISSION  AND  REFLECTION  FOR  CRYSTALS  OF  N  LAYERS. 

Here  we  consider  the  reflection  and  transmission  of  x  rays  through  a  crystal  composed  of 
N  layers.  The  calculation  uses  MDP  to  add  all  multiple  reflections  which  occur  in  the 
crystal.  The  solutions  are  calculated  in  terms  of  the  transmission  per  layer  x. 
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Wavelength(A) 
(c)  X  =  6.053  A 


Figure  C-2.  Quartz  1010  rocking  curves  under  Bragg  reflection 
for  both  the  n  and  <r  components. 


x  =  (-1)mexp(-ri), 


(C.20) 


where 

r|  =  ±yjss  -(g  +  e)2  .  (C.21) 

The  parameters  s,  cr  and  s  are  defined  as  in  the  case  of  a  reflection  from  a  thick  crystal. 
The  sign  in  ri  is  chosen  such  that  the  real  part  of  r\  is  negative  so  the  wave  is  attenuated  as 
is  passes  through  the  crystal. 

The  reflected  intensity  and  transmitted  intensities  for  jV  layers  are  given  by 

>|2 

(C.22) 


Rn  - 


$on 

to 

1 

*0 

1 

<N 

K 

1 

1 _ 

To 

To 

_1-(S0  /T0)2x2AJ 

and 


lN 


Ton 

2 

xN 

To 

X- 

1-(S0/T0)2 

1-(S0/T0)2x2N 


(C.23) 


C.2.1  Transmission  and  Reflection  from  N  Layers  of  Quartz  1010. 

The  transmission  and  reflection  for  quartz  10K)  was  calculated  for  N  =  1000,  N  =  5000, 
and  TV  =  10000  under  the  Bragg  condition  X  =  2d  sinQ.  The  calculations  were  performed 
for  both  states  of  polarization  at  X  =  5.739  A.  The  values  of  MF(0)  and  MF(B)  for  Quartz 
1010  at  this  wavelength  are  MF(0)  =  2.86  x  1016  +  4.786  x  1015  i  and  MF(Q)  =  -5.80  x 
1015-  1.320  x  1015z. 
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The  results  of  the  transmission  and  reflection  calculations  for  the  different  number  of 
layers  are  shown  in  Figure  C-3.  Also  included  in  these  plots  is  the  reflected  intensity  for 
the  semi-infinite  crystal  for  comparison.  The  plots  show  that  away  from  the  Bragg  angle, 
the  transmission  through  the  crystal  is  the  same  as  that  from  bulk  material  absorption  of  x 
rays.  This  is  true  for  both  states  of  polarization  showing  that  the  x  rays  that  pass  through 
the  crystal  which  are  not  at  the  Bragg  angle  are  not  polarized. 

In  the  vicinity  of  the  Bragg  angle,  the  reflection  for  the  a  polarization  increases  with 
crystal  thickness  as  expected  until  it  reaches  the  semi-infinite  crystal  results.  Here  the 
transmission  has  a  decrease  from  the  bulk  absorption  transmission  value.  This  is 
reasonable  because  as  more  intensity  is  reflected,  less  should  be  transmitted.  On  the  other 
hand,  the  n  polarization  shows  very  little  reflection.  This  is  occurs  because  the  Bragg 
angle  is  near  45°  where  the  n  polarization  is  not  reflected  with  any  significant  intensity. 
The  corresponding  transmission  shows  little  variation  from  the  bulk  transmission  near  the 
Bragg  angle. 
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Wavelength(A) 


Figure  C-3.  Transmission  Tn  and  reflection  Rn  of  5.739  A  from  Quartz  1(H0 
for  N  layers  and  reflection  R  from  a  semi-infinite  crystal. 


APPENDIX  D 
RAY  TRACING 


D.l.  RAY  TRACING  FOR  JOHANN  SPECTROMETERS. 

Ray  tracing  involves  tracing  each  x  ray  from  the  source  to  the  detector  to  get  an  intensity 
distribution  at  the  detector.  In  the  calculations  for  the  Phoenix  experiment,  we  were 
interested  in  performing  ray  tracing  for  spectrometers  in  the  Johann  configuration  for 
both  a  point  and  line  source. 

The  Johann  spectrometer  consists  of  a  crystal  bent  to  radius  of  curvature  of  R  with  a 
detector  placed  on  a  circle  of  diameter  R.  A  point  or  line  source  is  placed  a  distance  D 
from  the  center  of  the  crystal  at  the  Bragg  angle  corresponding  to  the  line  being  observed. 
Figure  D-l  shows  this  setup. 

The  procedure  for  finding  the  resulting  intensity  at  the  detector  is  as  follows.  Find  the 
flux  at  the  detector  at  a  particular  wavelength  for  a  point  in  the  source.  Sum  the 
contributions  of  flux  from  all  points  in  the  source  at  this  given  wavelength.  This  was  the 
calculation  performed  by  Morita  to  obtain  the  resolving  power  of  the  spectrometer  for  a 
line  source  at  a  given  wavelength/31^  We  extended  this  calculation  by  summing  over  the 
spectral  profile  at  each  point  to  produce  a  complete  spectral  profile  from  the  source.  The 
ray  tracing  equations  derived  from  Morita  are  presented  as  the  starting  point. 
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Figure  D-l.  Johann  Spectrograph. 
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We  assume  that  we  have  a  source  with  a  normalized  spectral  profile  J(k  -  Xq,x ,y)  and 
radiance  H(x,  y).  For  a  small  area  in  the  source  ds  dz,  the  flux  emitted  into  the  solid 
angle  d't  d\y  in  the  directions  of  the  ray  being  traced  is  given  by 

J{X-kQ, x,y)H(x, y)  cos(^)  dX  ds  dz dt,  dy ,  (D.l) 

where  %  is  the  angle  relative  to  the  normal  of  ds  in  the  direction  of  the  crystal  as  shown 
in  Figure  D-2  and  dy  is  the  angle  of  elevation  of  the  ray  out  of  the  plane  of  the 
spectrometer.  The  total  flux  received  at  the  detector  for  a  point  source  at  a  single 
wavelength  is  then  given  by 

F  =  /(A,)|  W(Q  -QB)cos(^)dXdsdzd^d\y  ,  (D-2) 

where  W(6-QB)  is  a  weight  corresponding  to  the  rocking  curve  of  the  crystal.  This  is 
the  reflectivity  of  the  crystal  as  a  function  of  crystal  angle  for  a  particular  wavelength. 
For  most  flat  crystals  the  rocking  curves  have  been  measured  to  be  Gaussian.  In  this 
calculation  we  assume  a  Gaussian  rocking  curve.  This  will  not  be  a  bad  assumption  since 
the  radius  of  curvature  of  the  crystal  used  in  the  experiment  is  large.  The  normalized 
rocking  curve  is  expressed 

=  (d.3> 

with 


e  *  Ihl, 

Qh  V  2tt 


(D.4) 
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where  Bfo  is  the  halfwidth  of  the  Gaussian  line  in  radians.  0  h  can  be  calculated  using 
Bragg’s  law  from  the  known  resolving  power  of  the  crystal  XI  AX. 

Using  the  relationship  for  arc  length  LSC(jd\sf  =  dh ,  with  Lscci  the  source  to  crystal 

distance  plus  crystal  to  detector  distance,  we  can  write  the  flux  in  terms  of  the  height  dh 
at  the  detector 

F=Cdzdhl(X)^^— ^-cos(^)dsd^.  (D.5) 

L scd 

For  a  point  source  s  with  a  wavelength  X ,  the  flux  .S'  at  the  detector  will  then  be 

Sxs  00  Js  (^)J  cos(£)  d\ ,  (D.6) 

Lscd 

with  Is  (X)  being  the  relative  intensity  of  the  input  point  at  the  given  wavelength  X. 
Numerically,  the  integral  in  ^  becomes  a  sum  over  discrete  angular  steps  starting  from 
the  right  end  of  the  crystal  and  proceeding  to  the  left  end  of  the  crystal.  Each  angle  is 
traced  and  the  flux  arriving  at  the  detector  is  added  to  the  appropriate  bin  in  the  detector. 
The  angular  step  size  must  be  small  enough  so  that  the  continuous  nature  of  the 
integration  is  not  lost.  Using  the  approximation  that  cos£,  «  1,  which  is  valid  because  the 
distance  to  the  crystal  is  large  compared  to  the  size  of  the  crystal,  the  discrete  form  of 
Eq.  (D.6)  becomes 

Ss(X)Kl,(Xj£-£—.  (D.7) 

£  t'scd 
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To  calculate  the  spectra  N(X)  for  an  extended  source  at  a  single  wavelength,  we  sum 
over  all  points  in  the  source 

^  ^  _  WAX) 

AWcc£,s,(X)  =  £4(r£-ii— .  (D.8) 

5-  s  E,  ^ scd 

The  final  profile  for  the  source  is  computed  by  summing  over  all  the  wavelengths  in  all  of 
the  source  points.  For  a  spectral  profile  JS(K  -  A,0) ,  the  flux  at  the  detector  is 


W 

*  0 ) «  £  mm = £  £  7,  p.  -  x  0  >2  .  (D.9) 

X  X  s  E,  L'scd 

This  sum  represents  the  total  flux  from  every  point  and  wavelength  in  the  source. 

To  generate  a  profile  we  must  add  the  flux  from  each  ray  in  the  source  to  the  appropriate 
place  on  the  detector.  To  do  this  we  have  to  determine  where  each  ray  strikes  the 
detector.  The  calculation  begins  by  finding  the  point  where  each  ray  strikes  the  crystal, 

Xc  and  7C  as  a  function  the  angle  <f>  from  the  x  axis  to  the  ray  being  traced  and  the  source 
point  Xs  and  Ys.  Referring  to  Figure  D-2  we  can  calculate^  from 


Xe  =z~xs  ~ J cos(O) - Xs , 


(D.10) 


by  using 


and 


tan((|)  R  +  A4>)  =  tan(O)  =  — 

Xc 


CD-11) 


R2  =j2  +g2 -2jgcos(W  +  0),  (D.12) 

where 
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g2  =X2+(R-YS)\ 


SP  =  cos  1  ( 


Jxf+(R-Y2)2 


), 


(D.13) 


and 


cos^F  +  O)  =  cos('F)  cos(O)  -  sin('F)  sin(O) 


-  Xs  cos(®)  +  (R  ~  Ys )  sin(O) 
"  Xs2+(R-Ys)2 


Solving  the  quadratic  Eq.  (D.12)  for  j 
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j  =  Xs  cos(d>)  +  (R-YS) sin(®) 

±cos(0)V(^  005(0)  +  ^- rjsin(O))2 -(X2  +(R-y)2 -R2).  (D.15) 

Plugging  this  expression  into  Eq.  (D.10)  we  then  have 

Xc  =  cos2($)[(l  +sec2  0)XS  +(R- FJtan(0)] 

±cos2(0)yjx2  +2Xs(R-Ys)tm®  +  (R-Ys)2  tan2  O  +  sec2  0(2i?F,  -X2  -Y2)  .  (D.16) 
With  the  appropriate  trigonometric  identities  this  can  be  rewritten 
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(D-17) 


Xc  =  cos2  <I>[(75  -  R  -  Xs  tan  ®)  tan  O 


By  defining  A  =  -  tan(<l>)  and  B  =  YS-  AXS  -  R , 

(D.18) 

From  Figure  D-3,  since  R2  =  P2  +  X2  and  Yc  =  R  -  P ,  we  find  Yc  =  R  -  ^R2  -X2  . 

Once  we  have  found  out  where  each  ray  strikes  the  crystal,  we  can  proceed  to  find  the 
position,  Xr  and  Yn  where  the  rays  strike  the  detector.  We  can  find  Xr  and  Yr  by  finding 
the  intersection  of  the  line  of  the  reflected  ray  from  the  crystal  with  the  Rowland  circle. 
From  Figure  D-2  we  see  that  line  from  where  the  ray  strikes  the  crystal  to  the  detector  is 
given  by  Y  =  (X  -  Xc )  tan  a  +  Yc .  Plugging  this  expression  for  Y  into  the  equation  of  the 
Rowland  circle,  given  by  X2  +  [Y  -  (R  /  2 )]2  =  (R  /  2)2 ,  with 

( Yr-  R  /  2)2  =  ( Xr  tana  +  Yc  -Xc  tana  -  R  /  2)2 

=  (Xr -2Xc)Xr  tan2  a  +  (2 YcXr  -  RXr -2 XCYC  +  RXC) tana2 

D 2 

-2XcXr  tan2  a  +  Yc2  -  RYC  +  — ,  (D.  19) 

we  have 

(1  +  tan2  a)X2  +  (27c  -  2XC  tana  -  R)Xrtma  +  Y2  -  RYC  +  RXC  tana 


-AB±^A2B2  -(1  +  A2)(B2 -R2) 


Figure  D-3.  Ray  intersection  with  perpendicular  and 
Rowland  circle  detectors. 


-2XC  Yc  tan  a  +  X2  tan2  a  =  0 .  (D.20) 

Simplifying  this, 

X2  sec2  a+(2 Yc  -2XC  tana -R)Xr  tana + if  -RYC  +RXC  tma-2XcYc  tana 

+X2  tan2  a  =  0 ,  (D.21) 

and  letting 

C  =  (2  Yc  -  2XC  tan  a  -  R)  tan  a  (D.22) 

and 

D  =  (Xc  tan  a  -  2YC  +  R)XC  tan  a  +  (Yc  -  R)YC ,  (D.23) 

we  then  have 
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(D.26) 


Bin  Number  =  f^er  BinsJ^  +  2£j  _ 

where  Z,  is  the  length  of  the  detector.  With  this  expression,  when  F  is  L  /  2,  the  flux  is 

added  to  the  last  bin  corresponding  to  the  number  of  bins  and  when  F  is  -L  /  2,  the  bin 

number  is  zero.  The  minus  sign  refers  to  spectra  on  the  detector  closest  to  the  crystal. 

Using  this  method  the  number  of  bins  can  be  changed  easily  if  necessary.  The  number  of 

bins  used  in  the  simulations  presented  here  is  100.  Using  too  many  bins  in  the 

calculation  can  present  problems.  The  main  problem  is  that  intermediate  bins  may  not 

get  filled  up.  This  is  caused  by  too  large  an  angle  increment  across  the  crystal  compared 

to  the  distance  between  bins.  This  can  be  resolved  by  decreasing  the  angular  steps  in  the 

ray  tracing  calculation,  but  this  results  in  increased  computational  time. 

With  the  rays  properly  assigned  to  their  bins,  then  it  is  a  simple  matter  finding  the 

distance  x  of  these  bins  on  the  detector 

-  Number  -  {Number  of  Bins  /  2)] 

Number  of  Bins 

When  the  bin  number  is  zero,  x  will  be  equal  to  -L  /  2  and  when  the  bin  number  equals 
the  number  of  bins,  x  will  be  equal  to  L  /  2. 

The  bins  then  have  to  be  converted  from  distance  on  the  detector  to  wavelength.  Figure 
D-4  shows  the  setup  required  for  the  calculation  of  the  detector  on  the  Rowland  circle. 
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circle  detector. 
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To  find  the  wavelength  corresponding  to  a  distance  F  from  the  center  of  the  detector  we 
must  find  the  new  Bragg  angle  corresponding  to  the  position  on  the  film. 

From  Figure  D-4  the  new  Bragg  angle  is 

=^-A  =  Qb+j  =  sin-I(“)  + 1 ,  (D.28) 

and  the  wavelength  corresponding  to  a  distance  F  on  the  Rowland  circle  detector  is  then 
simply  given  by 

=  2d  sin(sin_1  {—)  +  .  (D.29) 

2  a  R 

Using  this  expression,  replacing  F  by  x,  we  can  calculate  the  wavelength  of  each  of  the 
bins  that  the  rays  are  assigned  to. 
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APPENDIX  E 

X-RAY  FILM  CALIBRATION 


E.l  CALIBRATION  OF  KODAK  BIOMAX  MR  FILM. 

Kodak  Biomax  MR  4099  film  was  used  in  the  spectrometers  for  the  Phoenix 
experiments.  This  film  is  a  high  speed,  high  contrast  film,  with  low  fog  designed  for  use 
in  DNA  sequencing  autoradiography.  The  film  is  coated  with  a  thin  gelatin  layer  to  help 
prevent  damage  to  the  emulsion.  Since  this  film  is  not  standardly  used  as  x  ray  film  for 
the  wavelength  region  below  10  A,  calibration  was  necessary  to  determine  the  saturation 
level  and  the  range  of  optical  densities  where  the  response  of  this  film  is  linear.  This 
calibration  was  performed  for  the  Lp  line  of  A1 XIII  using  a  vacuum  spark  source. 

The  calibration  was  done  by  controlling  the  exposure  of  the  film  and  measuring  the 
resulting  optical  density  of  the  film.  The  technique  of  Troussel  and  Dumont  as  shown  in 
Figure  E-l  was  employed  for  this  task.(40)  A  flat  crystal  spectrograph  was  used  to 
disperse  the  x  rays  from  an  x  ray  source  source  and  a  stack  of  varying  thickness 
absorption  foils  were  placed  in  front  of  the  film  to  vary  the  exposure  across  the  spectral 
lines. 

The  source  of  the  x  rays  used  for  the  calibration  of  the  Biomax  MR  film  was  a  vacuum 
spark  with  an  aluminum  anode  and  cathode.  Since  the  yield  of  x  rays  on  vacuum  spark  is 
low,  several  hundred  shots  had  to  be  taken  to  produce  the  intensities  required  to  calibrate 
the  film. 
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X-Ray  Be  Filter  Crystal  Wedge  ;  Film 

Source  Y 

X 


Figure  E-L  Expenmentai  arrangement  for  film  calibration  (Dumont  and  Troussei). 
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A  flat  crystal  spectrograph  utilizing  a  single  ADP  101  crystal  with  2d  =  10.648  A  was 
employed.  The  spectrograph  varied  slightly  in  design  from  the  spectrograph  presented  in 
appendix  A  in  that  only  one  crystal  was  used.  The  spectrograph  was  designed  in  such  a 
way  that  the  ray  reflected  from  the  center  of  the  crystal  would  strike  the  center  of  the  film 
perpendicular  to  it.  The  crystal  was  located  a  distance  of  32.4  cm  from  the  source.  The 
distance  from  center  of  the  film  to  the  crystal  was  6  cm  and  the  angle  the  crystal  made 
with  the  center  of  the  source  was  38.9°,  in  order  to  place  the  Lp  line  in  the  center  of  the 
film.  This  configuration  resulted  in  a  dispersion  of  2.2  A/mm  on  the  film.  The  aluminum 
spectra  with  no  filters  obtained  with  the  flat  crystal  spectrograph  for  the  vacuum  spark  is 
shown  in  Figure  E-2.  Table  E-l  contains  the  line  assignments  for  the  spectra. 


T able  E- 1 .  Aluminum  spectra  from  vacuum  spark. 


Line  ID 

State 

Transition 

Wavelength  (A) 

Designation 

1 

ALXIII 

ls-7p 

5.493 

Lyman  \ 

2 

ALXIII 

ls-6p 

5.534 

Lyman  e 

3 

ALXIII 

ls-5p 

5.605 

Lyman  8 

4 

ALXIII 

ls-4p 

5.739 

Lyman  y 

> 

AL  XIII 

ls-3p 

6.053 

Lyman  p 

6 

AL  XII 

1  s2- 1  s5p 

6.1745 

NA 

For  the  calibration  to  work  correctly  a  scan  across  the  Lp  line  needed  to  have  a  flat  region 
in  intensity,  so  that  the  difference  in  intensities  getting  to  the  film  was  only  due  to 
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5.0  5.2  5.4  5.6  5.8  6.0  6.2  6.4 

Wavelength  (A) 

(a)  Aluminum  spectra  from  vacuum  spark 


Distance  (mm) 

(b)  Density  scan  across  Lp  line  (5) 


Figure  E-2.  Aluminum  spectra  obtained  from  vacuum  spark  and  density 
scan  across  the  Lp  line  (5). 
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the  absorption  foils  and  not  to  a  variation  in  source  intensity.  A  scan  across  the  Lp  line  is 
shown  in  Figure  E-2.  For  the  region  designated,  the  intensity  does  not  vary  more  than 
10%.  The  filters  were  then  placed  so  the  steps  were  within  this  region.  The  absorption  foils 
placed  in  front  of  the  film  consisted  of  a  stack  of  six  1-mil  Be  foils  offset  such  that  the 
width  of  each  step  was  0. 1  mm,  resulting  in  foil  thicknesses  which  ranged  from  1  to 
6-mil. 

The  film  was  scanned  using  a  Joyce  Loebl  scanning  densitometer  (Appendix  F).  An 
aperture  of  0. 1  was  used  with  an  occult  of  0.7  mm,  and  a  scan  width  of  50  pm.  Either  the 
C383  wedge  or  H5Z  wedge  was  used  depending  upon  the  maximum  density  of  the  film. 

The  density  scale  on  the  densitometer  was  calibrated  on  each  film  scan  by  using  at  least 
two  of  the  following  filters:  Melles  Griot  03FNG085  ( D  =  0.5),  03FNG015  (£>  =  1.0),  and 
03FNG023  ( D  =  2.0).  To  find  the  saturation  level  of  the  film,  a  piece  of  film  was 
overexposed  and  scanned.  This  resulted  in  a  saturation  density  of  2.73  ±  0.03. 

The  scans  of  the  Lp  lines  where  the  step  filters  were  used  are  shown  in  Figure  E-3.  To 
analyze  the  data,  the  background  fog  level  (typically  0.39  ±  0.02)  was  subtracted.  The 
recorded  density  of  the  various  steps  was  tabulated  for  each  set  of  data.  The  standard  x 
ray  transmission  calculations  for  various  thicknesses  of  Be  were  used  to  scale  the 
theoretical  exposure  of  the  film.(27)  The  simple  empirical  model  of  Dumont  et.  al.  was 
then  used  to  generate  response  curves.  The  expression  from  Dumont  et  al.  for  the  optical 
density  is  given  by 
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Figure  E-3.  Scans  across  the  L  p  lines  with  Be  step  filters. 
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where  L  is  an  arbitrary  exposure  defined  where  an  exposure  of  L0  leads  to  an  optical 
density  of  Dsa/2.  L0  was  chosen  to  be  500.  The  equation  for  optical  density  was  then  fit 
to  all  of  the  data  sets.  For  each  set  a  k  value  was  obtained.  The  data  with  the  fits  are 
shown  in  Figure  E-4.  The  value  of  k  varied  from  1.33  to  1.58.  The  discrepancy  of  it 
between  data  sets  can  be  attributed  to  poor  filter  placement  and  the  difficulty  of 
determining  the  intensity  at  each  filter  step.  Poor  filter  placement  resulted  in  a  variation 
of  the  intensity  across  the  film  in  addition  to  the  variation  imposed  by  the  filters.  The 
problem  of  determining  the  intensity  levels  at  each  step  resulted  from  two  factors,  the 
quality  of  the  step  filter  and  the  noise  caused  by  graininess  in  the  film.  The  poor  quality 
of  the  step  filter  came  from  the  difficulty  of  producing  the  Be  foil  stack  without  chipping 
away  the  edges  of  the  brittle  Be  foils.  A  higher  quality  step  filter  could  reduce  this 
problem  and  produce  better  intensity  plateaus.  The  noise  could  be  reduced  by  enlarging 
the  distance  between  steps  so  that  the  densitometer  can  scan  with  a  wider  slit.  When 
using  a  vacuum  spark,  however,  this  makes  it  more  difficult  to  get  all  of  the  steps  in  the 
region  of  constant  intensity  in  the  scan  across  the  line. 

From  the  data  it  can  be  seen  that  the  region  of  linear  response  is  approximately  1.7  orders 
of  magnitude.  It  is  necessary  to  make  sure  the  lines  analyzed  in  the  experiment  do  not 
have  peak  optical  densities,  with  fog  subtracted,  that  exceed  this  level. 
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Figure  E-4.  Film  calibration  curves  for  L  p. 


APPENDIX  F 
DATA  SCANNING 


F.l  J O Y CE-LOEBL  DENSITOMETER. 

The  spectral  lines  were  scanned  from  the  film  using  a  Joyce  Loebl  microdensitometer, 
which  measures  the  opacity  of  the  film  as  a  function  of  position.  The  schematic  diagram 
for  the  densitometer  is  given  in  Figure  F-l .  Light  from  an  incandescent  bulb  (A)  is  split 
into  two  optical  paths.  The  light  from  one  path  passes  through  the  sample  (J)  and  the 
light  from  the  other  path  is  sent  through  a  calibrated  wedge  with  linearly  varying  optical 
density  (K).  The  beams  are  brought  together  at  a  photomultiplier  tube  (E)  which  is 
switched  to  look  alternately  at  each  beam.  A  servo  motor  (N)  moves  the  optical  wedge  so 
that  the  signal  from  both  beam  legs  are  the  same.  The  recording  pen  is  coupled  to  the 
wedge  to  produce  the  displacement  in  the  pen. 

F.1.1  Densitometer  Parameters  for  Phoenix  Data  Scans. 

The  alignment  procedure  for  the  sample  follows  from  the  instruction  manual  for  the 
Joyce  Loebl  instrument.  For  the  Phoenix  data,  a  5x  objective  with  0.1  aperture  lens  was 
used  in  the  instrument.  The  iris  was  set  so  that  the  slit  height  on  the  sample  was  0.7 
mm.  A  source  slit  width  of  25  was  used  which  corresponds  to  50  pm  on  the  sample.  An 
arm  ratio  of  2  was  used  for  all  Johann  Spectrometer  data  and  5  for  all  flat  crystal  data. 

The  scans  were  run  with  a  forward  speed  of  3.5.  The  H2Z  wedge  was  used  for  the 
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Figure  F-l.  Joyce  Loebi  scanning  densitometer. 
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analysis  of  the  Johann  spectrometer  data  and  the  C383  wedge  was  used  for  the  flat  crystal 
data. 


F.1.2  Data  Digitization. 

For  data  analysis  it  was  necessary  to  have  the  data  in  digital  form.  The  analog  data  from 
the  Joyce  Loebl  were  digitized  using  an  Apple  color  scanner  in  combination  with  the 
program  Datathief  1 .0.8,  which  extracts  the  data  from  the  scan.  The  scanning  procedure 
follows.  First  the  film  was  scanned  on  the  densitometer  in  the  usual  way.  A  reference 
line  was  placed  on  the  film  horizontal  to  the  spectra  with  a  fixed  length  of  200  mm.  The 
densitometer  spectra  was  then  scanned  using  Ofoto  2.0,  a  standard  scanning  program  for 
the  Apple  Macintosh  computer.  The  scanning  software  was  set  to  scan  a  grayscale  photo 
with  8:G  scan  bits,  135  scan  dpi,  8:G  print  bits,  and  135  print  dpi.  Once  the  spectra  had 
been  scanned  it  was  checked  for  alignment.  The  horizontal  reference  marker  was 
checked  to  ensure  that  it  was  aligned  with  the  page.  The  spectra  was  rescanned  if  the 
scan  was  not  aligned  correctly.  The  data  was  then  rotated  to  bring  the  spectra  in  the 
proper  orientation.  The  focusing  feature  of  Ofoto  2  was  employed  on  the  data  to  raise  the 
contrast  of  the  data  so  that  it  was  more  easily  read  by  the  DataThief  program.  Finally  the 
data  was  saved  in  the  PICT  ColorSync  Image  format  for  use  by  the  DataThief  software. 
The  final  step  to  digitizing  the  data  was  to  extract  the  numbers  from  the  scanned  image  of 
the  spectra.  The  program  DataThief  was  used  for  this  purpose.  The  program  reads  the  x 
and  y  components  of  a  line  into  numbers  in  text  format.  Once  an  image  has  been  loaded 
ir.ro  the  program  the  x,  y  scales  have  to  be  assigned.  For  the  Phoenix  data,  the  x  scale  was 
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set  using  the  200  mm  reference  marker.  This  calibrated  the  scale  of  the  spectra  to  mm 
(On  the  film  the  actual  scale  was  half  of  this  since  the  arm  ratio  used  was  2).  They  scale 
was  chosen  at  the  extremes  of  the  page.  A  few  parameters  had  to  be  set  before  the  data 
was  finally  read  in.  The  trace  width  and  trace  height  tell  the  program  how  many  pixels  to 
look  horizontally  or  vertically  if  there  is  a  gap  in  the  data,  since  the  data  from  the  scanner 
will  not  be  perfect.  A  trace  width  of  4  pixels  and  trace  height  of  30  pixels  was  used. 

With  the  auto  trace  feature  turned  on,  the  spectral  lines  were  scanned  and  data  saved  as 
text  files.  The  data  was  then  ready  for  manipulation  by  a  standard  scientific  program 
package. 
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